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Abstract.  The  distribution  of  age-of-air  in  a  two-dimensional,  zonally  averaged 
model  of  the  atmosphere  is  examined.  The  age-of-air  is  determined  from  mixing 
ratios  of  a  fictitious  tracer  with  a  linearly  increasing  tropospheric  source.  The 
impact  of  horizontal  mixing  is  isolated  by  performing  parallel  calculations  with 
and  without  horizontal  mixing  terms  in  the  tracer  evolution  equation.  We  examine 
the  evolution  of  the  tracer  over  several  model  years,  using  different  tropical  heat 
sources  as  well  as  different  amounts  of  planetary  wave  mixing  and  drag.  Typical 
ages  in  the  midstratosphere  when  horizontal  mixing  is  included  range  from  3  to 
over  5  years,  depending  primarily  on  the  strength  of  the  model’s  prescribed  tropical 
heat  source.  These  values  of  age  are  within  the  range  of  observational  estimates 
based  on  long-lived  tracers  with  increasing  tropospheric  mixing  ratios,  such  as  HF, 
CO2  and  SF6.  Age  determined  from  the  horizontally  unmixed  tracer  is  reduced 
by  1  to  2  years.  However,  the  evolution  of  the  unmixed  tracer  is  consistent  with 
observed  behavior  of  tropical  water  vapor  oscillations  in  data  from  the  HAlogen 
Occult  at  ion  Experiment  (HALOE).  We  found  stronger  than  expected  sensitivity 
of  middle  atmospheric  age  to  prescribed  heating  in  the  tropical  troposphere.  This 
sensitivity  in  our  model  arises  from  small  resulting  changes  to  the  residual  vertical 
velocity  in  the  tropics.  These  vertical  velocity  changes  are  in  turn  associated  with 
net  heating  changes  of  0.05-0.10  K  d"1  in  the  tropical  lower  stratosphere.  We  argue 
that  changes  of  this  magnitude  are  within  the  uncertainties  imposed  by  natural 
variations  in  the  tropics. 


1.  Introduction 

Accurate  stratospheric  measurements  of  long-lived 
tracers  with  increasing  tropospheric  mixing  ratios  have 
recently  become  available.  These  include  SFg  [Elkins 
et  al.  1996;  Harnisch  et  al,  1996],  CO2  [Schmidt  and 
Khedim,  1991;  Boering  et  al ,  1994],  and  HF  [Russell 
et  al,  1996].  Mixing  ratios  of  such  tropospheric  source 
gases  can  be  used  to  estimate  the  amount  of  time  re¬ 
quired  for  air  parcels  to  travel  from  the  troposphere 
to  other  locations  in  the  atmosphere.  This  is  done 
straightforwardly  by  comparing  the  mixing  ratio  of  a 
tropospheric  source  gas  at  say  40  km  altitude  with  a 
time  history  of  its  tropospheric  mixing  ratio.  The  time 
delay  between  the  occurrence  of  a  given  mixing  ratio  in 
the  troposphere  and  its  appearance  at  higher  altitudes 
is  referred  to  as  the  “age”  of  the  air  mass,  or  simply 
age  of  air.  This  concept  was  first  introduced  by  Bischof 
et  al  [1985]  to  interpret  differences  in  middle  strato¬ 
spheric  and  tropospheric  mixing  ratios  of  CO2.  Values 
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of  age  determined  in  this  way  have  been  found  to  vary 
between  3  and  7  years  in  the  midstratosphere  (25-35 
km)  [Schmidt  and  Khedim ,  1991],  Hall  and  Prather 
[1993]  and  more  recently  Waugh  et  al  [1997]  exam¬ 
ined  the  distribution  of  age  in  the  stratosphere  in  three- 
dimensional  general  circulation  model  (GCM)  simula¬ 
tions.  Hall  and  Prather  found  ages  between  3  and  4 
years,  while  Waugh  et  al.  found  much  older  ages,  from 
6  years  to  8-10  years  depending  on  the  parameteriza¬ 
tion  of  gravity  waves  which  was  employed.  These  stud¬ 
ies  used  a  single  year’s  output  from  the  GCMs  to  per¬ 
form  multiyear  tracer  simulations,  so  sensitivity  of  age 
to  variations  in  dynamic  or  thermodynamic  inputs  to 
the  GCMs  were  not  assessed  in  detail.  However,  it  is 
clearly  of  interest  that  such  large  differences  in  age  exist 
between  atmospheric  models.  Hall  and  Plumb  [1994]  ex¬ 
amined  the  effects  of  the  time  behavior  of  tropospheric 
sources  on  determinations  of  stratospheric  age. 

Age-of-air  in  the  middle  atmosphere  has  also  been 
defined  as  the  time  needed  by  a  Lagrangian  parcel  to 
travel  from  the  equatorial  troposphere  to  another  loca¬ 
tion  in  the  middle  atmosphere  [e.g.  Rosenlof ,  1995]. 
Using  parcels  advected  by  a  2-D  residual  circulation 
inferred  from  United  Kingdom  Meteorological  Office 
(UKMO)  analyses  of  temperature,  Rosenlof  obtained 
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ages  of  less  than  3  years  in  most  of  the  stratosphere  from 
Lagrangian  travel  times.  Older  values  were  obtained  in 
the  polar  vortices.  However,  rapid  horizontal  motions  in 
the  atmosphere  may  produce  significant  differences  be¬ 
tween  the  Lagrangian  age  derived  from  2-D  circulations 
and  the  mean  age  of  an  ensemble  of  parcels  in  a  fully  3-D 
flow  .  This  reflects  the  fact  that  an  air  parcel  in  the  2-D 
circulation  is  actually  composed  of  subelements  which 
can  have  widely  diverging  reverse  trajectories  back  to 
origin  at  the  surface  or  injection  into  the  stratosphere. 
Hall  and  Plumb  [1994]  introduced  the  concept  of  an  “age 
spectrum”  to  encompass  the  variety  of  paths  taken  by 
air  parcels  between  the  troposphere  and  any  particular 
point  in  the  middle  atmosphere.  The  mean-age,  deter¬ 
mined  from  long-lived  tracers  with  linearly-increasing, 
tropospheric  mixing  ratios,  is  equal  to  the  first  moment 
of  Hall  and  Plumb’s  age  spectrum  [Hall  and  Plumb, 
1994].  However,  irreversible  horizontal  mixing  broad¬ 
ens  the  age  spectrum  so  that  its  first  moment  is  not 
equal  to  2-D  Lagrangian  travel  times.  Thus  we  should 
expect  age  of  air  based  on  mixing  ratio  meaurements  of 
long-lived  trace  gases  to  be  different  than  Lagrangian 
ages  based  on  the  2-D  residual  circulation. 

The  rate  of  transport  from  the  troposphere  into  the 
middle  atmosphere  is  determined  by  a  balance  between 
residual  mean  upward  advection  in  the  tropics  and  hori¬ 
zontal  exchange  of  air  between  tropics  and  midlatitudes 
[e.g.,  Boering  et  al,  1995;  Plumb ,  1996;  Volk  et  al , 
1996].  Comparing  age  of  air  based  on  mixing  ratio  with 
age  determined  from  Lagrangian  parcel  calculations  will 
shed  light  on  the  realism  of  this  balance  in  2-D  models. 
In  zonally  averaged  2-D  models,  irreversible  horizontal 
mixing  by  atmospheric  motions  is  typically  represented 
as  a  diffusion  process  [e.g.,  Plumb  and  Ko ,  1992;  Garcia 
et  al ,  1992],  The  accuracy  of  horizontal  mixing  par am- 
eterizations  is  a  major  and  well-recognized  uncertainty 
in  2-D  models. 

On  the  other  hand,  the  calculation  of  radiative  net 
heating  is  generally  regarded  as  a  strong  point  of  zon¬ 
ally  averaged  models,  since  it  can  be  accurately  calcu¬ 
lated  using  a  number  of  radiatively  active  gases.  There¬ 
fore  “tuning”  of  less  well  understood  drag  parameters  in 
these  models  can,  in  principle,  be  accomplished  by  ob¬ 
taining  accurate  simulations  of  wind  and  temperatures 
in  the  middle  atmosphere  in  the  presence  of  known  dis¬ 
tributions  of  radiatively  active  gases.  When  both  cor¬ 
rect  drag  and  heating  distributions  are  present  in  the 
model,  a  correct  residual  mean  meridional  circulation 
is  guaranteed.  This  procedure  is  described  by  Summers 
et  al  [1997]  and  Siskind  et  al  [1997]  for  the  current 
Naval  Research  Laboratory  (NRL)  2-D  model. 

Unfortunately,  in  the  troposphere,  nonradiative  pro¬ 
cesses  such  as  sensible  surface  heating  and  latent  heat¬ 
ing  by  water  vapor  condensation  are  the  dominant  in¬ 
puts  of  heat.  These  processes  are  not  generally  well  rep¬ 
resented  in  zonally  averaged  2-D  models.  These  uncer¬ 
tainties  would  be  of  little  interest  to  middle  atmosphere 
modelers  if  their  effects  were  restricted  to  the  tropo¬ 
sphere  and  lower  stratosphere.  However,  in  practice, 
we  have  found  that  model  simulations  of  N2O  are  sen¬ 
sitive  to  the  strength  of  the  specified  heat  source  in  the 


troposphere  [Summers  et  al ,  1997].  Net  heating  rates 
in  the  tropical  lower  stratosphere  are  particularly  sensi¬ 
tive  to  variations  in  infrared  radiation  welling  up  from 
below  [e.g..  Eluszkiewicz  et  al,  1996].  Also,  it  is  now 
thought  that  vertical  velocity  at  the  tropical  tropopause 
and  in  the  lower  stratosphere  may  be  sensitive  to  or  even 
controlled  by  lower  stratospheric  drag  in  midlatitudes 
[e.g.,  Iwasaki ,  1992;  Rosenlof  and  Holton ,  1993;  Yulaeva 
et  al ,  1994;  Holton  et  al,  1995].  Thus,  vertical  velocity 
in  the  tropical  lower  stratosphere  may  be  sensitive  to 
numerous  processes  which  are  poorly  or  incompletely 
incorporated  in  2-D  models  and  as  will  be  seen  in  this 
study,  ascent  rates  through  the  lower  stratosphere  ex¬ 
ert  an  important  control  on  the  age-of-air  distribution 
in  the  middle  atmosphere. 

In  this  study  we  will  present  equilibrium  distributions 
of  age  of  air  from  a  2-D  zonally  averaged  model  of  the 
atmosphere.  We  will  use  the  NRL  chemical/dynamical 
model  [Bacmeister  et  al ,  1995;  Summers  et  al,  1997] 
to  calculate  the  evolution  of  a  fictitious  long-lived  tracer 
with  a  linearly  increasing  tropospheric  mixing  ratio. 
The  mixing  ratio  of  this  tracer  in  the  middle  atmosphere 
will  be  used  to  obtain  an  equivalent  age  which  should 
correspond  to  those  obtained  from  long-lived  trace  gases 
such  as  SF6  or  CO2.  We  will  perform  parallel  age  cal¬ 
culations  using  a  tracer  that  is  immune  to  horizontal 
mixing.  Age  from  this  “unmixed”  tracer  should  re¬ 
semble  that  which  would  be  obtained  from  Lagrangian 
parcel  calculations  using  the  model  residual  circulation. 
There  are  two  results  from  this  study:  (1)  we  demon¬ 
strate  that  2-D  models  with  parameterized  planetary 
wave  drag  and  mixing  can  obtain  an  approximately  cor¬ 
rect  balance  between  upward  advection  in  the  tropics 
and  horizontal  mixing,  and  (2)  we  show  that  variations 
in  prescribed  heating  of  the  tropical  troposphere  can 
lead  to  1  to  2  year  differences  in  the  age  of  the  model 
middle  atmosphere.  Point  2  is  principally  a  modeling 
issue,  since  the  tropospheric  dynamics  of  a  zonally  aver¬ 
aged  model  of  this  kind  are  not  realistic.  Nevertheless, 
in  exploring  the  sensitivity  of  the  model  to  tropospheric 
heating  we  found  that  much  of  the  middle  atmosphere 
response  to  tropospheric  heating  variations  was  asso¬ 
ciated  with  small  variations  in  the  net  heating  of  the 
tropical  lower  stratosphere.  These  variations  in  strato¬ 
spheric  net  heating  rates  are  within  the  uncertainties 
present  in  current  radiative  transfer  calculations,  due 
to,  for  example,  the  aerosol  loading  of  the  stratosphere. 

The  outline  of  this  paper  is  as  follows:  In  section  2  we 
briefly  describe  the  formulation  of  the  NRL  2-D  model. 
We  also  describe  the  current  state  of  model  parame- 
terizations  for  gravity  wave  drag,  radiative  transfer  and 
other  physical  processes.  In  Section  3  we  discuss  results 
from  several  simulations.  After  presenting  our  base¬ 
line  model  climatology  and  age  distribution  in  Sections 
3.1  and  3.2,  we  will  examine  the  impact  of  changing 
the  planetary  wave  forcing  (section  3.3)  as  well  as  the 
prescribed  tropospheric  heating  (section  3.4)  on  model 
age-of-air.  In  section  4  we  examine  in  more  detail  the 
response  of  the  model  meridional  stream  function  to 
tropospheric  heating  parameters.  Finally,  in  section  5 
we  summarize  our  results. 
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2.  Description  of  Model  and  Numerical 
Experiments 

A  basic  description  of  the  dynamical  portion  of  the 
NRL  2-D  model  is  given  by  Bacmeister  et  al.  [1995].  A 
brief  summary  of  the  model  formulation  is  given  here. 
The  model  is  based  on  the  Transformed  Eulerian  Mean 
(TEM)  formulation  [e.g.,  Andrews  et  al.,  1987].  Evolu¬ 
tion  equations  for  both  zonally  averaged  potential  tem¬ 
perature  0  and  angular  momentum  M  are  solved: 


3M 

dt 


■  +  v' 


3M  *  3M 
ad4>+W  ~8z 


=  X 


(la) 


86  86  ,86 

8i  +  v  ~^  +  w  Yz~n 


(lb) 


where  zonally  averaged  angular  momentum  density  is 
given  by 

M  =  Ua  cos(0)  +  fia2  cos2(0)  (lc) 


The  numerical  advection  of  M  and  9  in  (la)  and  (lb) 
is  accomplished  using  the  scheme  described  by  Prather 
[1986].  This  scheme  is  highly  nondiffusive  and  nondis- 
persive  and  compares  well  with  analytical  solutions  of 
the  advection  equation  [Prather,  1986;  Rood ,  1987;  Shia 
et  al ,  1990]  The  term  X  denotes  the  total  momen¬ 
tum  forcing  from  all  sources,  including  planetary  wave 
Eliassen-Palm  flux  (EP-flux)  divergence  and  gravity 
wave  drag.  The  term  H  denotes  the  total  thermody¬ 
namic  forcing  from  all  sources.  The  individual  com¬ 
ponents  of  X  and  H  are  discussed  below.  All  other 
symbols  have  the  same  meaning  as  in  the  work  of  An¬ 
drews  et  al.  [1987]  unless  otherwise  noted.  The  ther¬ 
mal  wind  relation  is  used  with  (la)  and  (lb)  to  form  a 
diagnostic  elliptic  equation  for  the  residual  meridional 
stream  function  0*.  The  stream  function  is  used  to 
calculate  the  meridional  and  vertical  residual  veloci¬ 
ties  v *  and  w* .  The  integration  of  (1)  and  the  solu¬ 
tion  of  the  meridional  stream  function  equation  includ¬ 
ing  a  geostrophic  adjustment  procedure  is  described  by 
Bacmeister  et  al.  [1995].  A  major  modification  to  the 
model,  since  Bacmeister  et  al.  [1995],  has  been  the  ad¬ 
dition  of  nonzero  phase  speed  gravity  waves  to  the  grav¬ 
ity  wave  parameterization  scheme.  The  current  status 
of  this  and  other  parameterizations  is  described  briefly 
next. 


2.1.  Momentum  Forcing 

The  total  momentum  forcing  X  in  the  model  can 
be  broken  down  into  components  representing  various 
physical  processes, 

3  3 

X  —  DgW  “  ®RayU  +  ~Q^{Xmorn  — {/)  +  V  •  Fpw  (2) 

We  now  briefly  describe  each  term  in  (2)  and  its  current 
implementation  in  the  model. 

2.1.1.  Gravity  wave  drag.  The  term  Dqw  rep¬ 
resents  drag  produced  by  dissipating  gravity  waves.  It 
is  the  only  component  of  X  that  has  been  substantially 
modified  since  Bacmeister  et  al.  [1995].  Drag  and  mix¬ 
ing  from  gravity  waves  with  non-zero  phase  speeds  have 


been  added  to  obtain  improved  simulations  of  the  meso¬ 
sphere  and  upper  stratosphere.  To  calculate  drag  and 
mixing  caused  by  gravity  waves  with  non-zero  phase 
speeds,  the  parameterization  uses  a  simplified,  hydro¬ 
static  wave  packet  approach  [Bacmeister,  1993;  Sum¬ 
mers  et  al.,  1997].  In  the  current  simulations  we  spec¬ 
ify  time  and  latitiude  independent  surface  amplitudes 
and  phase  speeds  for  a  six  wave  spectrum  of  traveling 
gravity  waves.  These  are  given  in  Table  1.  For  drag 
from  zero  phase  speed  gravity  waves  the  simple  cubic 
law  drag  of  Bacmeister  et  al.  [1995]  is  retained.  The 
scheme  is  described  in  more  detail  by  Summers  et  al. 
[1997].  _ 

2.1.2.  Rayleigh  drag.  The  second  term  £*RayJ7 
is  an  ad  hoc  Rayleigh-damping  term.  This  is  a  crudely 
represented  sink  of  angular  momentum  in  the  model 
stratosphere  and  troposphere.  It  is  assumed  to  account 
for  missing  or  imperfectly  parameterized  dynamics  in 
the  model,  such  as  low-altitude  orographic  drag  from 
large- amplitude  mountain  waves,  or  drag  from  high- 
wavenumber  planetary  waves.  In  the  present  study  we 
will  use  QfRay=0.01  d-1  for  all  calculations. 

2.1.3.  Vertical  mixing.  The  third  term 
3z(KmomdzU)  represents  vertical  redistribution  of  mo¬ 
mentum  by  turbulent  and  molecular  diffusion.  Its  ef¬ 
fects  on  our  simulations  are  negligible  below  the  upper 
mesosphere,  where  large  values  of  Km0m  can  be  associ¬ 
ated  with  gravity  wave  breaking  and  where  large  ver¬ 
tical  shears  in  zonal  winds  may  be  present  due  to  the 
effects  of  gravity  wave  drag. 

2.1.4.  Planetary  Wave  Drag.  The  term  V-Fpw 
represents  EP-flux  divergence  associated  with  dissipat¬ 
ing  or  breaking  planetary  waves.  Planetary  wave  effects 
are  parameterized  as  in  the  work  of  Garcia  [1991],  by 
forcing  a  time  dependent,  linear  planetary  wave  model 
with  zonal  winds  and  temperatures  from  the  zonally 
averaged  TEM  model.  In  this  study  we  specify  a  time- 
independent,  surface  geopotential  height  perturbation 
for  wavenumber  1  described  by, 


#(4>,zolt)=< 


0 

Znh  |  sin  30/2  | 
Zsh  |  sin  30/2  | 
0 


if  0  >  60° 
if  0°  <  0  <  60° 
if  -60°  <  0  <  0° 
if  0  <  -60° 


Table  1.  Parameters  for  nonzero  phase  speed  gravity 
waves  used  in  all  simulations  discussed  in  this  study. 
Column  1  gives  wave  phase  speed.  Column  2  gives 
peak  vertical  displacement  amplitude  associated  with 
the  waves  at  the  launch  altitude  of  %5.0  km.  Column 
3  gives  the  number  of  wave  packets  with  these  phase 
speeds  and  amplitudes  encountered  around  a  latitude 
band. 


Phase  Speed  ms1 

Amplitude  m 

Number  per  Band 

+20 

40.0 

5.0 

-20 

40.0 

5.0 

+40 

30.0 

3.0 

-40 

30.0 

3.0 

+40 

15.0 

3.0 

-40 

15.0 

3.0 
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Table  2.  Experiments  to  be  discussed  in  the  text.  Column  1  indicates  the  number 
of  the  experiment  in  the  text.  Column  2  gives  the  magnitude  of  the  prescribed 
tropospheric  heating  used  in  the  experiment.  Column  3  indicates  whether  inputs  to 
the  thermodynamic  forcing,  such  as  O3  or  NMC  surface  temperatures,  have  been 
“symmetrized”  with  respect  to  hemisphere  and  season.  Columns  4  and  5  indicate  the 
amplitude  of  the  planetary  wave  forcing  in  the  northern  and  southern  hemispheres. 
Experiment  1  represents  the  baseline  simulation.  Results  for  experiment  2  are  similar 
in  many  respects  to  those  for  experiment  1,  except  that  hemispheric  asymmetries  are 
supressed.  In  experiments  2,  3,  and  4,  symmetric  planetary  wave  forcing  amplitudes 
are  varied,  while  tropospheric  heating  magnitude  Ho  is  held  fixed.  In  experiments  5 
and  6,  asymmetric  planetary  wave  forcing  amplitudes  equal  to  those  in  the  baseline 
simulation  (experiment  1)  are  used,  while  the  magnitude  of  prescribed  tropospheric 
heating  Ho  is  varied.  In  experiment  7,  the  tropospheric  heating  is  the  same  as  in 
experiment  1,  but  a  small  ~0.3  K  d"1  heat  source  centered  at  25  km  is  added. 


Experiment 

Ho  K  d"J 

Inputs  to  H 

Zm  m 

Zsh  m 

1 

3.0 

Baseline  simulation 
asymmetric 

350 

175 

2 

3.0 

Vary  Planetary  Waves  With  Fixed  Ho 
symmetric 

250 

250 

3 

3.0 

symmetric 

125 

125 

4 

3.0 

symmetric 

350 

350 

5 

2.0 

Vary  Ho  With  Fixed  Planetary  Waves 
asymmetric 

350 

175 

6 

4.0 

asymmetric 

350 

175 

7 

3.0  +  H$tT&t 

Add  Stratospheric  Heating 
asymmetric 

350 

175 

We  retain  the  possibility  to  select  different  surface  forc¬ 
ing  amplitudes  in  the  southern  and  northern  hemi¬ 
spheres  as  a  crude  way  of  modeling  observed  hemi¬ 
spheric  asymmetries  in  planetary  wave  activity.  Ta¬ 
ble  2  gives  values  of  Zsh  and  Znh  and  other  param¬ 
eters  for  each  of  the  seven  experiments  discussed  in  this 
study.  Earlier  studies  have  used  Zsh&  0.5^nh  to  ob¬ 
tain  reasonable  hemispheric  asymmetries  in  such  fea¬ 
tures  as  wintertime  trace  gas  distributions  within  the 
polar  vortices  [ Garcia  et  a/.,  1992;  Bacmeister  et  a/., 
1995].  The  planetary  wave  EP-flux  divergence  is  also 
used  to  calculate  a  mixing  coefficient  Kpw  to  represent 
rapid  isentropic  horizontal  mixing  by  breaking  plane¬ 
tary  waves,  thought  to  be  the  most  important  cause  of 
rapid  horizontal  mixing  in  the  stratosphere  [e.g.,  McIn¬ 
tyre  and  Palmer ,  1983].  A  more  detailed  description  of 
the  planetary  wave  model  is  given  by  Bacmeister  et  al 

[1995]. 

2.2.  Thermodynamic  Forcing 

The  contributions  to  the  thermodynamic  forcing  H  in 
(lb)  can  also  be  broken  down  into  those  from  a  number 
of  different  physical  processes, 

W  =  Wgas  4*  ^latent  +  H&traX  “  abdy(^  —  $NCEp) 

+£(*“”4*)  (3) 

The  significance  of  each  term  and  its  current  parame¬ 
terization  is  described  below. 


2.2.1.  Radiative  transfer.  The  first  term  7fgae 
is  the  net  heating  due  to  exchange  of  solar  and  ter¬ 
restrial  radiation  by  gases  in  the  atmosphere.  We  use 
two  radiative  transfer  parameterizations  to  calculate  in¬ 
frared  heating  and  cooling  in  different  layers  of  the  at¬ 
mosphere.  Climatological  O3,  CO2,  and  H2O  are  used. 
Details  are  discussed  by  Summers  et  al  [1997]. 

2.2.2.  Prescribed  tropospheric  or  “latent” 
heating.  The  second  term  Tfiatent  is  an  ad  hoc  pre¬ 
scribed  heating  function.  It  is  intended  to  account  for 
nonradiative  thermal  forcing  by  condensation  and  rain- 
out  of  tropospheric  water  vapor  as  well  as  other  pro¬ 
cess  such  as  sensible  heat  transfer  from  the  surface. 
The  effects  of  changing  the  distribution  and  strength 
of  this  prescribed  heating  will  be  a  particular  focus  of 
this  study.  We  will  use  the  basic  distribution  shown  in 
Figure  1  for  the  prescribed  tropospheric  heat  source  or 
“latent  heating.”  A  similar  approach  to  modeling  the 
nonradiative  heating  of  the  troposphere  is  used  in  other 
2-D  models  [e.g.,  Garcia  et  al,  1992]. 

Our  prescribed  tropospheric  heating  in  Figure  la  is 
a  simple  analytical  function  of  latitude  and  height  de¬ 
scribed  by 


^latent (<£,  M)  =  'Ho  [Fi(<j>,Z,t)  +  -F2(<M)] 
where 


■^1  (^)  zi  0 

=  MAX  [l-(^fef)2(3^)2,0 


(4a) 


(4b) 
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Figure  1.  (a)  Tropospheric  heat  source  or  ‘latent  heating”  as  a  function  of  latitude  and 

height  in  October.  Heating  is  shown  for  Wo=3Kd’1,  the  value  used  in  most  of  the  experiments 
discussed.  In  experiments  5  and  6,  Hq  is  varied,  but  the  heating  distribution  in  Figure  la  is 
unchanged  except  for  proportional  changes  in  strength,  (b)  Additional  stratospheric  heat  source 
which  is  used  in  experiment  7. 


and, 


F^,z)  =  g(z)ex p  (- 


|0|  ~  0mid 
A^mid 


) 


(4c) 


Thus,  The  function  in  (4a)  is  the  sum  of  a  deep,  broad 
tropical  heat  source,  represented  by  F\ ,  and  a  shallower 
midlatitude  source,  represented  by  J^.  The  width  of  the 
tropical  source  in  degrees  latitude  is  given  by  A^trop- 
We  fix  A^trop  to  35°  for  all  the  calculations  discussed 
here.  The  center  of  the  tropical  heating  is  given  by 
^trop(O)  which  oscillates  between  10°S  on  June  21  and 
10°  N  on  December  21.  The  horizontal  scale  of  the  mid¬ 
latitude  heat  source  A0mid  is  set  to  30°.  The  center 
^mid  is  fixed  at  45°  in  both  hemispheres  for  the  entire 
year.  The  vertical  distribution  of  midlatitude  heating 
is  given  by  g(z)>  which  we  take  to  be  a  sharply  peaked 
function  with  a  maximum  value  of  1  at  5  km  and  linear 
fall  off  to  zero  within  2.5  km  on  both  sides  of  the  peak. 
The  depth  of  the  tropical  heat  source  is  given  by  dtrop 
and  is  set  to  15  km  for  all  of  experiments  in  this  study. 
The  strength  of  the  heat  source  7io  will  vary  within  the 
range  2.0  K  d*"1  to  4.0  K  d”1  (Table  2).  The  motiva¬ 
tion  for  the  prescribed  heating  distribution  in  Figure  1 
comes  from  inferred  climatologies  of  atmospheric  latent 


heating  [e.g.,  Newell  et  al.,  1970].  These  studies  show 
a  deep  tropical  source  of  latent  heat  as  well  as  shallow 
“wings”  in  midlatitudes.  Typical  maximum  values  of 
latent  heating  from  these  climatologies  lie  between  2.5 
Kd"1  and  3.5  K  d”1.  Characteristic  depths  for  tropical 
latent  heating  are  around  15  km. 

2.2.3.  Lower  stratospheric  tropical  heating. 
The  term  7fstrat  is  included  to  incorporate  small  geo¬ 
physical  and  theoretical  uncertainties  in  the  net  heating 
of  the  tropical  lower  stratosphere.  Geophysical  uncer¬ 
tainties  include  those  arising  from  aerosols  [Kinne  et  ah, 
1992;  Rosenfield  et  al,  1997],  or  cloud  albedo  variations 
[e.g.,  Eluszkiewicz  et  al.,  1996].  Theoretical  uncertain¬ 
ties  in  radiation  algorithms  arise  for  example  from  the 
suite  of  radiating  gases  considered  [e.g.,  Eluszkiewicz  et 
al,  1996].  In  all  of  the  simulations  to  be  discussed  here 
except  for  Experiment  7  (Table  2),  7f8trat  is  zero  every¬ 
where.  For  experiment  7  we  use 


Tfstrat 


=  0.3 


K  d^(Ad) 


where  <£troP  is  as  in  (4b),  A<£strat  is  set  to  25°,  d8trat  is 
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set  to  25  km,  and  Azstrat  is  set  to  10  km.  7fstrat=0  is 
imposed  outside  of  the  range  |^-<£trop(OI<A^trop,  and 
\z— dstrat|<Azstrat.  This  function  is  shown  in  Figure 
lb.  The  peak  of  7fatrat  occurs  immediately  above  that 
of  Wutent  at  all  times  and  thus  also  oscillates  between 
10°N  and  10°S. 

2.2.4.  Surface  boundary  heating  and  cooling. 

The  term  c*bdy(0  —  0ncep)  prevents  the  lower  tropo¬ 
sphere  from  cooling  to  excessively  cold  temperatures 
at  the  poles  and  in  winter  midlatitudes.  It  is  nonzero 
only  in  the  lowest  vertical  level  of  the  model.  At  this 
level,  the  model  potential  temperature  is  relaxed  to  ob¬ 
served  potential  temperatures  derived  from  a  seven  year 
average  of  National  Center  for  Environmental  Predic¬ 
tion/Climate  Analysis  Center  (NCEP/CAC)  analyses. 
This  term  can  be  thought  of  as  a  crude  parameterization 
of  surface  heating  effects,  including  poleward  transport 
of  heat  by  the  oceans. 

2.2.5.  Vertical  mixing.  The  last  term  in  (3)  rep¬ 
resents  vertical  redistribution  of  heat  by  turbulent  and 
molecular  diffusion.  Turbulent  diffusion  of  heat  is  a  sig¬ 
nificant  contribution  to  the  total  thermodynamic  forc¬ 
ing  in  both  the  tropical  troposphere  and  in  the  meso¬ 
sphere.  In  the  mesosphere,  high  values  of  ifthrm  arise 
from  breakdown  of  gravity  waves.  In  the  troposphere 
we  prescribe  large  values  for  Zfthrm  in  the  tropics  to 
coincide  approximately  with  the  tropical  latent  heat 
source.  Peak  values  of  ATthrm  in  the  tropical  troposphere 
are  taken  to  be  105  cm2  s-1  at  <t>=<j>trop  and  to  fall  off 
quadratically  in  both  height  and  latitude.  The  scales 
for  /<thrm  are  15  km  in  the  vertical  and  35°  in  latitude. 
These  large  vertical  diffusivitites  act  on  strong  vertical 
gradients  of  potential  temperature  and  result  in  signif¬ 
icant  fluxes  of  heat.  This  distribution  of  Ffthrm  will  be 
fixed  for  all  of  the  experiments  described  here. 


2.3.  Age  of  Air  and  Experimental  Design 


The  mean  age  of  air  (henceforth  denoted  by  T  follow¬ 
ing  Hall  and  Plumb  [1994])  in  our  model  simulations  will 
be  determined  by  tracking  the  evolution  of  a  fictitious 
“age  tracer”  pi{<j>,  z ,  t)  with  a  tropospheric  source  P  and 
no  losses  anywhere.  The  tracer  evolution  equation  is 
straightforward, 


d±  +  v*^-  +  w*dJL 

dt  ad<f>  dz 


=  P  + 


cos<t>Kyv^  +  PK» If  (5) 


As  was  the  case  for  (la)  and  (lb)  the  numerical  ad- 
vection  of  the  tracer  is  performed  using  the  scheme  of 
Prather  [1986].  The  tracer  (age  tracer)  is  initially  uni¬ 
form  throughout  the  atmosphere.  Its  initial  value  is  set 
to  100  units.  Below  10  km  we  set  P  to  be  a  constant, 
latitudinally  uniform  source  with  a  strength  of  7  units 
yr-1,  and  above  10  km  we  set  P= 0.  The  present  cal¬ 
culations  include  a  background  or  “floor”  value  for  the 
horizontal  eddy  diffusion  Kyy  of  109  cm2  s_1.  Eddy  dif¬ 
fusion  from  planetary  wave  breaking  ATpw  is  used  wher¬ 
ever  it  exceeds  the  background  value;  that  is,  Kyy  = 


MAX[(109  cm2  s-1),  ifpw]*  The  background  value 
for  Kyy  is  included  to  account  for  horizontal  mixing  in¬ 
duced  by  dynamical  processes  in  the  atmosphere  not  in¬ 
cluded  in  our  wave  1  planetary  wave  model.  These  could 
include  higher-order  planetary  waves,  quasi-geostrophic 
turbulence,  or  inertiogravity  waves  [e.g.,  Pierce  et  a/., 
1994].  The  vertical  mixing  coefficient  for  tracers  Kzz 
is  set  equal  to  Tfthrm ,  the  vertical  eddy  diffusivity  for 
heat. 

In  order  to  examine  the  effects  of  horizontal  mixing 
on  model  transport  we  also  track  the  evolution  of  a  re¬ 
lated  tracer  /i*.  This  tracer  is  initialized  and  forced 
identically  to  pi.  However,  when  (5)  is  used  to  calculate 
the  evolution  of  p*,  we  set  Kyy=0  identically.  Verti¬ 
cal  mixing  is  retained  in  the  evolution  of  p* .  Thus  by 
comparing  p  and  p*  the  maximum  possible  impact  of 
horizontal  mixing  on  tracer  evolution  in  the  model  can 
be  assessed. 

Transport  and  vertical  mixing  in  the  model  tropo¬ 
sphere  prevent  the  mean  tropospheric  mixing  ratio  of  p 
or  p*  from  possessing  exactly  the  behavior  specified  by 
our  choice  of  P.  We  integrate  the  model  for  12  years 
and  perform  a  linear  fit  to  the  equatorial  mixing  ra¬ 
tio  of  the  age  tracers  at  the  model  lower  boundary.  In 
practice,  this  low-level  mixing  ratio  is  a  nearly  linear 
function  of  time  with  a  slope  of  around  5.9  units  yr-1. 
The  result  of  the  linear  fit  is  then  used  to  calculate  the 
age  at  a  particular  time  and  location  according  to, 

p  _  £?  0  ~  £?  Q 

(A/i/A  *)m  W 

where  p  can  be  either  the  horizontally  mixed  tracer  or 
the  unmixed  tracer  p* .  Depending  on  whether  pi  or  //* 
is  used  in  (6)  we  refer  either  to  a  “mixed”  age  V  or  to 
an  “unmixed”  age  T*  for  the  model. 

The  mixed  age  T  determined  from  (6)  should  be 
equivalent  to  age  estimated  from  measurements  of  long- 
lived  tracers  with  time-dependent  tropospheric  sources, 
while  the  unmixed  age  T*  should  be  close  to  estimates 
of  age  based  on  air-parcel  travel  times  in  a  2-D  La- 
grangian  calculation  such  as  that  shown  by  Rosenlof 
[1995].  Some  differences  between  T*  and  Lagrangian 
travel  times  could  still  arise  due  to  vertical  mixing  and 
implicit  numerical  diffusion  that  may  still  affect  the  evo¬ 
lution  of  p* .  However,  for  most  of  the  stratosphere  Kzz 
is  small  and  numerical  diffusion  in  our  advection  scheme 
is  known  to  be  small  [  Prather ,  1986;  Shia  et  al ,  1991; 
Bacmeister  et  a/.,  1995].  Therefore,  we  expect  the  un¬ 
mixed  age  based  on  p*  to  be  reasonably  close  to  the  true 
2-D  Lagrangian  travel  times  that  would  be  obtained  us¬ 
ing  our  model’s  meridional  circulation.  In  a  2-D  TEM 
model,  T  and  T*  as  defined  in  (6)  can  be  expected  to 
be  quite  different  The  magnitude  of  difference  may  be 
related  to  the  width  of  the  age-spectrum  discussed  by 
Hall  and  Plumb  [1994]  and  Hall  and  Waugh  [1997]  since, 
at  least  in  the  tropics,  T*  and  the  the  time-of- arrival  for 
the  peak  in  the  age  spectrum  (the  “modal  time”  Hall 
and  Waugh  [1997])  should  be  similar. 
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ZONAL  WIND  AND  TEMPERATURE  DURING  SOLSTICE 
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Figure  2.  Zonal  wind  and  temperature  as  a  function  of  latitude  and  height  at  both  solstices 
for  experiment  1,  the  baseline  simulation.  The  top  panels  show  model  zonal  winds;  (left)  close  to 
northern  winter  solstice,  January  6,  and  (right)  close  to  summer  winter  solstice,  July  6.  Hatching 
and  dashed  contours  show  easterly  winds.  Thick  contours  are  drawn  every  20ms"1,  thin  contours 
every  10  m  s-1.  The  bottom  panels  show  temperatures;  (left)  near  northern  winter  and  (right) 
near  southern  winter  solstice.  Hatching  shows  temperatures  below  160  K.  Thick  contours  are 
drawn  every  40  K,  thin  contours  every  10  K. 


3.  Results 

3.1.  Model  Climate 

Figure  2  shows  zonal  wind  and  temperature  at  sol¬ 
stice  conditions  for  a  simulation  with  peak  troposo- 
pheric  heating  ft0=3  K  d-1  and  planetary  wave  forcing 
Zsh=  175  m  and  Znh=  350  m.  As  discussed  by  Sum¬ 
mers  et  al.  [1997],  upper  stratospheric  and  mesospheric 
zonal  wind  and  temperature  fields  look  reasonable  when 
compared  with  the  climatologies  of  Fleming  et  ai  [1988] 
and  Randel  [1992].  Minimum  summer  mesopause  tem¬ 
peratures  are  somewhat  lower  in  the  present  simulations 
than  in  those  discussed  by  Summers  et  al.,  Tm in~  120  K 
compared  with  Tm in~140  K.  Zonal  wind  reversals  in  the 
upper  mesosphere  are  somewhat  more  pronounced  in 
the  present  simulations.  These  relatively  minor  changes 
to  the  model  upper  mesosphere  are  due  to  differences  in 
the  settings  for  the  nonzero  phase  speed  gravity  waves 
shown  in  Table  1.  Summers  et  al.  use  a  hemispheri- 
cally  asymmetric  cubic  law  drag,  which  is  not  done  in 
the  present  simulations.  Here  latitudinally  uniform  cu¬ 
bic  law  drag  is  used.  This  results  in  somewhat  weaker 


zonal  wind  and  temperature  asymmetries  above  50  km 
than  were  the  presented  by  Summers  et  al  [1997]. 

In  the  lower  stratosphere  and  troposphere,  larger 
disagreements  between  model  and  observations  occur. 
Subtropical  jets  in  both  hemispheres  are  positioned  5° 
to  10°  too  close  to  the  equator  compared  with  NCEP 
analyses  [Randel,  1992]  and  Geophysical  Fluid  Dynam¬ 
ics  Laboratory  (GFDL)  analyses  [Peixoto  and  Oort, 
1992],  Observed  wintertime  subtropical  jets  in  both 
hemispheres  exhibit  maximum  winds  between  30°  and 
40°  latitude,  whereas  the  model  jets  peak  between  20° 
and  30°.  Peak  wind  speeds  in  the  winter  subtropical 
jets  are  observed  to  be  near  35  ms-1.  One  of  the  worst 
aspects  of  the  model  simulations  is  an  excessively  cold 
midlatitude  and  polar  tropopause.  Wintertime  extrat- 
ropical  tropopause  temperatures  in  the  model  are  10  to 
20  K  too  cold,  while  summertime  tropopause  tempera¬ 
tures  are  20  to  30  K  too  cold.  It  is  worth  noting  that 
in  simulations  without  planetary  wave  drag,  wintertime 
tropopause  temperatures  are  even  10  K  colder  than  in 
Figure  2.  This  suggests  that  our  large  errors  in  polar 
and  midlatitude  tropopause  temperatures  may  be  due 
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Amplitude  and  Phase  of  Annual  cycle  in  Temp. 


Latitude 


T  at  18.62  km 


Figure  3.  (top)  Amplitude  (solid  contours)  and  phase  (dashed  contours)  of  annual  cycle  in 
temperature  as  a  function  of  latitude  and  height  for  experiment  1.  Peak-to-peak  amplitudes  are 
twice  what  is  displayed,  (bottom)  Time  series  of  model  temperature  oscillation  at  18  km  at  2.5°S 
and  2.5°N.  Annual  average  has  been  subtracted. 


in  part  to  the  omission  of  heat  transport,  or  forced  de¬ 
scent,  by  baroclinic  eddies  or  other  dissipating  waves  in 
the  troposphere  and  lower  stratosphere.  Transport  of 
angular  momentum  by  baroclinic  eddies  may  also  move 
the  subtropical  jets  poleward 

3.1.1.  Tropical  tropopause  temperatures 
While  the  model  simulation  of  extratropical  tropopause 
temperatures  is  somewhat  deficient,  the  simulation  of 
tropical  tropopause  temperatures  is  quite  good.  Val¬ 
ues  vary  annually  between  196  K  and  200  K,  in  rea¬ 
sonable  agreement  with  observed  tropical  tropopause 
temperatures.  Figure  3  shows  amplitude  and  phase  for 
the  annual  component  of  temperature  variability  as  a 
function  of  latitude  and  height  for  the  present  exper¬ 
iment.  Model  tropopause  temperatures  in  the  tropics 
exhibit  an  annual  cycle,  with  minimum  values  during 
northern  hemisphere  (NH)  winter  and  maximum  values 
in  southern  hemisphere  (SH)  winter.  Figure  3  shows 
the  Fourier  amplitude  of  the  annual  cycle  in  the  model 
temperatures,  so  that  peak  to  peak  variations  are  twice 
the  value  displayed.  The  amplitude  of  the  Fourier  an¬ 
alyzed  annual  cycle  is  somewhat  weak,  2  to  3  K  peak 
to  peak,  compared  the  observed  4  to  10  K  range  re¬ 
ported  by  numerous  investigators  [e.g.,  Rosenlof  and 
Holton ,  1993;  Yulaeva  et  al,  1994].  However,  an  exam¬ 
ination  of  the  model  temperature  time  history  at  the 
tropical  tropopause  (Figure  3,  bottom)  shows  that  the 


total  difference  between  minimum  and  maximum  trop¬ 
ical  tropopause  temperatures  is  close  to  4  K. 

The  phase  of  the  model’s  annual  cycle  of  tempera¬ 
ture  at  the  tropical  tropopause  is  in  excellent  agree¬ 
ment  with  observations  cited  by  Yulaeva  et  al.  with 
temperatures  peaking  between  August  1  and  September 
1  at  18  km.  However,  at  lower  altitudes  the  observed 
phase  varies  rapidly  with  observed  tropical  tempera¬ 
tures  peaking  between  February  1  and  March  1  below 
15  km.  In  contrast,  the  phase  of  our  model  annual  cycle 
remains  relatively  constant  below  15  km.  This  diagree¬ 
ment  may  reflect  our  model’s  incomplete  simulation  of 
tropospheric  processes. 

Above  20  km  the  model’s  annual  cycle  in  tropical 
tropopause  temperatures  decreases  rapidly.  This  is  con¬ 
sistent  with  observations  [e.g.,  Yulaeva  et  al.,  1994; 
Garcia  et  al,  1997].  By  30  km,  variability  in  model 
tropical  temperatures  is  dominated  by  a  semiannual  os¬ 
cillation  (SAO).  This  is  also  in  accord  with  observations 
of  tropical  temperatures  [e.g.,  Garcia  et  al ,  1997].  Our 
model’s  SAO  is  discussed  in  greater  detail  by  Summers 
et  al  [1997]. 

3.2.  Basic  Distribution  of  T  and  T* 

We  now  describe  the  model  distribution  of  T  and  T* 
in  the  baseline  simulation.  A  latitude-height  section  of 
T  after  10.5  years  of  model  integration  is  shown  in  Fig- 
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Mixed 


Latitude 

age  jxn-ext  +  apw  Date=3l/l2 
age*  jxn-ext  +  apw  Date  =  31/12 

Figure  4.  Age  of  air  as  a  function  of  latitude  and 
height  on  day  31  of  model  year  12  (after  11.5  year  inte¬ 
gration).  (a)  Mixed  age  T,  (b)  unmixed  age  T*  . 


ure  4a.  The  corresponding  unmixed  age  F*  is  shown  in 
Figure  4b.  We  show  the  distribution  on  day  31  of  model 
year  11,  i.e.  northern  hemisphere  winter.  Both  F  and 
T*  have  a  characteristic  distribution  determined  largely 
by  the  Brewer-Dobson  (residual  mean  meridional)  cir¬ 
culation  of  the  atmosphere.  In  the  model  troposphere, 
prescribed  rapid  vertical  mixing  as  well  as  strong  ver¬ 
tical  motion  in  the  tropics  maintain  low  T  values  (<  1 
year).  However,  significant  latitudinal  gradients  in  T* 
develop  throughout  the  troposphere.  Higher  values  of 
T*  are  found  near  the  poles,  probably  as  a  result  of 
continual  injection  from  above  by  annual  mean  descent 
over  the  poles.  Nevertheless,  in  the  tropical  troposphere 
both  T  and  T*  have  values  below  1  year. 

In  the  tropical  stratosphere  between  15  and  40  km 
a  strong  vertical  gradient  is  present  in  both  mixed  and 
unmixed  ages.  This  gradient  is  present  primarily  be¬ 
cause  of  the  weak  mean  vertical  and  meridional  motion 
in  this  part  of  the  atmosphere.  Comparing  T  (Figure 
4a)  with  T*  (Figure  4b),  horizontal  mixing  actually  in¬ 
tensifies  this  vertical  gradient.  This  is  most  easily  inter¬ 
preted  as  a  simple  consequence  of  exchanging  descend¬ 
ing  parcels  with  low  values  of  fi  from  the  polar  and 
midlatitude  lower  stratosphere  with  ascending  tropical 
air  parcels  containing  higher  values  of  fi.  This  exchange 
leads  to  net  aging  of  the  tropical  stratosphere  and  all  of 
the  mesosphere,  compared  to  the  age  which  would  be 


present  in  the  absence  of  horizontal  mixing.  The  only 
part  of  the  atmosphere  where  mixing  leads  to  reduced 
age  is  the  extratropical  lower  stratosphere.  Above  40 
km,  strong,  annually  varying  meridional  winds  in  the 
tropics  along  with  strong  vertical  winds  at  the  poles 
rapidly  distort  and  homogenize  both  the  /i  and  the  /4* 
tracer.  This  leads  to  a  relatively  uniform  age-of-air 
distribution  everywhere  above  40  km.  The  mixed  age 
above  40  km  is  between  3.5  and  4.5  years,  while  the 
unmixed  age  is  between  2.0  and  2.5  years. 

Thus  the  basic  distribution  of  T  in  Figure  4,  and  to 
a  lesser  degree  T*,  consists  of  two  relatively  uniform 
regions,  one  in  the  troposphere  and  another  above  the 
midstratosphere.  These  two  regions  are  divided  by  a 
strong  vertical  gradient  in  age  in  the  lower  stratosphere, 
which  is  intensified  by  horizontal  mixing.  Further  ef¬ 
fects  of  planetary  wave  mixing  are  evident  in  the  step¬ 
like  structure  of  the  contours  of  T  in  the  lower  strato¬ 
sphere.  Going  south  from  90°N  at  around  20  km,  the  T 
contours  are  at  first  relatively  flat  with  a  value  around  4 
years.  Between  45°  and  55° N  a  relatively  sharp  drop  in 
T  to  values  less  than  3  years  is  encountered.  From  45° 
to  around  20° N  T  is  relatively  flat  varying  between  2.5 
to  2  years.  Between  15°  and  20° N  the  tropical  upward 
bulge  in  the  age  contours  begins.  At  the  equator,  ages 
of  1  year  or  less  are  found  up  to  25  km.  The  location 
of  the  well  mixed  region  between  15°  and  45° N  corre¬ 
sponds  roughly  to  the  location  of  the  midlatitude  surf 
zone,  which  is  bounded  by  the  polar  vortex  edge  and 
subtropical  barrier.  The  edge  of  the  subtropical  barrier 
is  somewhat  too  close  to  the  equator,  an  error  which  is 
probably  related  to  incorrect  position  of  the  subtropical 
jets  (Figure  2). 

The  general  features  of  the  lower-stratospheric  age 
distribution  shown  in  Figure  4a  are  in  good  agreement 
with  SF6  measurements  taken  from  the  ER-2  [Elkins 
et  a/.,  1996].  However,  values  of  T  at  35  km  at  mid¬ 
latitudes  are  low  compared  to  age  estimates  based  on 
balloon-borne  CO2  measurements  taken  between  1976 
and  1990  [Schmidt  and  Khedim ,  1991]  as  well  as  re¬ 
cent  balloon-borne  SF$  measurements  [Harnisch  et  a/., 
1996].  The  CO2  measurements  give  values  of  age  in  the 
midlatitude  midstratosphere  that  vary  between  2  and 

7  years  with  a  mean  of  5.6  years,  with  the  more  re¬ 
cent  CO 2  measurements  suggesting  a  somewhat  lower 
mean.  The  SFg  measurements  of  Harnisch  et  al.  yield 
very  high  ages  at  midlatitudes,  in  some  cases  exceeding 

8  years.  Our  baseline  T  is  also  somewhat  lower  than 
recently  published  estimates  of  age  based  on  HALOE 
HF  measurements  [Russell  et  al.  1996].  Russell  et  al. 
report  ages  5.9±2  years  at  50  km,  so  our  model  result 
of  ~4.0  years  is  low  but  still  within  Russell  et  al.’s  esti¬ 
mated  errors. 

The  unmixed  age  T*  is  not  an  observable  quantity. 
However,  T*  in  Figure  4b  is  in  rough  agreement  with  the 
age  shown  by  Rosenlof  [1995,  Figure  17]  except  in  the 
polar  lower  stratosphere  below  20  km  where  our  maxi¬ 
mum  F*  of  5  years  is  1.0  to  1.5  years  greater  than  that 
obtained  by  Rosenlof.  However,  above  25  km  and  in  the 
tropical  lower  stratosphere  our  T*  is  in  good  agreement 
with  the  results  of  Rosenlof.  The  age  shown  by  Rosenlof 
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is  based  on  travel  times  of  Lagrangian  air  parcels  ad- 
vected  by  a  residual  circulation  inferred  from  UKMO 
analyses  of  temperature  for  1992  and  1993.  It  is  of  inter¬ 
est  that  the  model  again  fails  most  visibly  in  the  polar 
lower  stratosphere.  Both  our  cold  bias  and  excessively 
high  unmixed  age  are  consistent  with  insufficient  dy¬ 
namically  forced  descent  in  the  lowest  part  of  the  polar 
stratosphere.  In  general  though,  the  similarity  between 
model  r  and  Rosenlof ’s  Lagrangian  age  calculations  is 
an  encouraging  sign  that  the  model  meridional  circula¬ 
tion  is  in  reasonable  agreement  with  that  inferred  from 
atmospheric  observations. 

3.3.  Sensitivity  to  Planetary  Wave  Forcing 

Planetary  wave  forcing  can  affect  the  model  age  dis¬ 
tribution  in  two  ways.  Planetary  wave  mixing  increases 
the  fraction  of  older  midlatitude  air  which  enters  the 
tropical  ascent  zone  of  the  model.  This  will  tend  to 
increase  the  age  of  the  middle  atmosphere.  However, 
increased  planetary  wave  momentum  forcing  can  also 
strengthen  the  meridional  circulation,  leading  to  more 
rapid  tropical  ascent  and  a  younger  middle  atmosphere. 
Unfortunately,  the  planetary  wave  forcing  in  this  model 
can  only  be  incompletely  specified  from  atmospheric  ob¬ 
servations.  Estimates  of  amplitudes  and  drag,  or  EP- 
flux  divergence,  for  stationary  planetary  wave  1  have 
been  made  from  data  [e.g,  Randel,  1992].  These  esti¬ 
mates  are  useful  as  checks  against  grossly  unrealistic 
planetary  wave  amplitudes.  However,  they  are  not  nec¬ 
essarily  useful  in  choosing  parameter  settings  that  yield 
good  simulations  of  zonally  averaged  winds  and  tem¬ 
peratures.  For  example  in  the  southern  winter  strato¬ 
sphere,  stationary  wave  1  is  observed  to  have  geopo¬ 
tential  amplitude  near  240  m,  which  is  far  below  50% 
of  the  northern  hemisphere  observed  amplitude  of  over 
960  m  [Randel,  1992].  Nevertheless,  we  obtain  rea¬ 
sonable  simulations  using  Zsn=0.5Z?m.  If  we  forced 
the  planetary  wave  model  with  something  like  the  ob¬ 
served  ratio  of  southern  to  northern  hemispheric  wave 
1  amplitudes  (~1 :4) ,  our  simulation  in  one  of  the  two 
hemsipheres  would  suffer.  We  speculate  that  this  is 
because  the  total  planetary  wave  effect  in  both  hemi¬ 
spheres  includes  higher  wavenumber  waves  and  nonsta¬ 
tionary  waves  which  do  not  possess  so  large  a  north- 
south  asymmetry  as  stationary  wave  1.  Our  model  cur¬ 
rently  simulates  only  planetary  wave  1  forced  by  orog¬ 
raphy.  The  effects  of  this  single  wave  are  expected  to 
account  for  those  of  a  complete  spectrum  of  stationary 
as  well  as  transient  planetary  waves.  It  is  therefore  of 
great  interest  to  examine  the  sensitivity  of  our  simula¬ 
tions  to  the  choice  of  planetary  wave  forcing  amplitudes 
Z$n  and  Znh-  In  order  to  restrict  parameter  space  and 
to  simplify  interpretation  of  the  results  in  this  Section, 
we  limit  our  attention  to  experiments  with  symmetric 
planetary  wave  forcing;  that  is,  Zsh=^nh=^nsh- 

Figure  5  shows  the  final  distributions  of  T  for  ex¬ 
periments  with  Znsh=125,  250,  and  350  m  (experi¬ 
ments  2-4,  Table  2).  The  final  T  for  experiment  2  with 
Znsh=250  m  (Figure  5a)  is  very  similar  to  that  for  ex¬ 
periment  1,  with  an  average  T  in  the  middle  atmosphere 


a  =  age*  jxn-ext+spw  Date  =  31/12 
b=age*  jxn-ext+spw_125  Date=31/12 
c  =  age*  jxn--ext+spw_350  Date  =  31/12 


Figure  5.  T  as  a  function  of  latitude  and  height 
on  day  31  of  model  year  12  for  experiments  with  fixed 
tropospheric  heating  and  varying  symmetric  planetary 
wave  forcing  amplitudes,  (a)  Znsh  =250  (experiment 
2).  (b)  Zksh= 125  (experiment  3).  (c)  Znsh  = 350  (ex¬ 
periment  4). 


close  to  4  years  in  both  experiments.  Generally  speak¬ 
ing,  the  broad  features  of  the  distribution  of  T  are  sim¬ 
ilar  in  all  of  the  experiments.  However,  T  in  the  strato¬ 
sphere  and  mesosphere  shows  a  systematic  dependence 
on  the  amplitude  of  planetary  wave  forcing.  Typical 
values  for  T  at  50  km  vary  from  a  low  value  of  3.5  years 
for  experiment  3  with  Znsh=125  m  (Figure  5b),  to  a 
high  value  4.5  years  for  experiment  4  with  Znsh=350 
m  (Figure  5c).  Thus,  there  is  a  systematic  decrease  in 
T  as  planetary  wave-forcing  amplitudes  are  reduced. 

In  contrast  to  the  systematic  variation  of  the  mixed 
age  T,  the  unmixed  stratospheric  age  T*  for  experiments 
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a  — age  jxn-ext  +  spw  Date  =  3l/12 
b  =  age  jxn-ext+spw_125  Date=31/12 
c  =  age  jxn-ext+spw_350  Date  =  31/12 

Figure  6.  Same  as  Figure  5  except  for  F\ 


2-4  in  Figure  6  reveals  little  consistent  variation  with 
Znsh-  The  average  value  of  T*  at  50  km  is  around 
2.5  years  for  all  of  the  experiments.  Minor  differences 
occur  between  experiments,  with  some  hint  of  reduced 
upward  motion  in  the  tropics  for  Znsh=350  m.  How¬ 
ever,  overall,  the  distribution  of  unmixed  age  in  all  three 
experiments  is  nearly  the  same.  This  can  only  be  the 
case  if  the  tropical  ascent  rates  in  all  three  experiments 
with  different  Znsh  are  similar.  Thus  the  increase  in  T 
with  increasing  planetary  wave  forcing  must  be  due  pri¬ 
marily  to  enhanced  horizontal  mixing.  Figure  7  shows 
time- averaged  Kyy  in  the  lower  stratosphere.  Zones  of 
strong  planetary  wave  mixing  in  midlatitudes  are  clearly 
evident,  and  the  magnitude  of  this  midlatitude  mixing 
depends  strongly  on  Znsh-  Little  if  any  enhancement 
of  Kyy  over  the  background  value  of  109  cm2  s_1  oc¬ 
curs  for  experiment  3  (dotted  line),  whereas  maximum 


Kyy  values  of  nearly  1010  cm2  s  1  occur  for  experiment 
4.  (dashed  line).  Figure  8  shows  profiles  of  time-  and 
latitude- averaged  residual  vertical  velocity  <  w **  >20 
in  the  tropics  for  experiments  with  varying  planetary 

wave  parameters.  In  the  following,  <  (  )  >4,  will  be 
used  to  indicate  a  time  average  over  3  model  years  and 
a  latitude  average  about  the  equator  between  +<j>  and 
—<j)  degrees.  Below  18  km,  <ttX>2o  has  a  clear  de¬ 
pendence  on  ^nsh  wifh  weaker  planetary  wave  forc¬ 
ing  giving  weaker  <uX> 20-  However,  little  effect  on 
<w*t  >20  is  noticeable  between  18  and  40  km,  where 
all  of  the  curves  lie  close  to  one  another. 

The  insensitivity  of  w*  to  increased  planetary  wave 
forcing  is  puzzling,  particularly  since  Kyy  and  EP-flux 
divergence  are  directly  connected  in  the  model’s  plane¬ 
tary  wave  parameterization.  Time  averaged  profiles  of 
V  •  Fpw  at  30°  N  are  shown  in  Figure  8.  As  Znsh  in¬ 
creases,  the  largest  increases  in  V  *  Fpw  are  seen  below 
22  km.  This  is  high  enough  to  create  a  mixing  zone  that 
can  interact  with  tropical  ascent.  However,  the  circula¬ 
tion  effects  of  this  added  drag  are  directed  downward, 
and  thus  have  little  impact  on  the  speed  of  tropical  as¬ 
cent  above  20-25  km.  Between  22  and  40  km,  V  •  Fpw 
at  30°N  for  Znsh=250  and  Znsh=350,  are  similar  in¬ 
dicating  that  planetary  waves  have  reached  saturation 
by  the  time  they  reach  these  altitudes.  As  pointed  out 
by  Garcia  [1991],  increased  planetary  wave  forcing  may 
be  reflected  primarily  in  the  size  of  the  planetary  wave 
breaking  region  rather  than  in  the  strength  of  the  drag 
in  a  model  with  saturation-limited  waves. 

This  pattern  of  response  to  changes  in  planetary  wave 
forcing  may  be  exaggerated  by  the  model’s  somewhat 
unrealistic  wind  distribution  in  the  troposphere.  The 
separation  of  tropospheric  and  middle  atmospheric  jets 
(Figure  2)  may  lead  to  early  saturation  of  planetary 
waves  originating  in  the  troposphere.  If  this  separa¬ 
tion  were  reduced,  changes  in  EP-flux  divergence  could 
be  distributed  over  a  deeper  layer  of  the  atmosphere. 
However,  winter  profiles  of  V  ■  Fpw  from  our  simula¬ 
tions  (Figure  8)  agree  reasonably  well  with  climatologies 
for  northern  hemisphere  winter  [e.g.,  Randel ,  1992]  al¬ 
though  our  maximum  winter  values  in  the  lower  strato¬ 
sphere  are  somewhat  lower  than  in  Randel,  which  could 
lower  the  sensitivity  of  uX  to  changes  in  planetary  wave 
forcing.  Nevertheless,  the  overall  wintertime  pattern  of 
EP  flux  in  our  simulations  agree  reasonably  well  with 
climatology.  In  section  4  we  explicitly  calculate  the  ver¬ 
tical  motion  forced  solely  by  mechanical  drag,  and  find 
it  to  be  small  compared  to  that  forced  by  other  terms. 
This  suggests  that  even  doubling  the  EP  flux  divergence 
in  our  simulations  would  have  a  limited  effect  on  tropi¬ 
cal  upwelling. 

3.4.  Sensitivity  to  Latent  Heating 

As  discussed  in  section  2  the  term  Wiatent  is  a  crude 
attempt  to  represent  the  heating  of  the  atmosphere  by 
irreversible  condensation  of  water  vapor.  Latent  heat¬ 
ing  is  as  yet  impossible  to  measure  directly,  so  it  must 
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Figure  7.  Multiyear  average  of  horizontal  eddy-mixing  coefficient  Kyy  between  15  and  22 
km  as  a  function  of  latitude  for  experiments  2-6.  Thick  solid  line  shows  Kyy  for  experiment  2. 
Dotted  line  shows  Kyy  for  experiment  3,  and  dashed  line  shows  curve  for  experiment  4.  Little 
enhancement  over  the  background  value  of  109cm2  s-1  occurs  for  experiment  3  where  Znsh^ 25 
m  is  used.  All  other  experiments  show  evidence  of  large  midlatitude  “surf  zones.”  Thin  lines 
connecting  symbols  show  curves  for  experiments  5  and  6,  i.e.,  fixed  asymmetric  planetary  wave 
forcing  with  varying  Wo*  Only  comparatively  small  changes  in  Kyy  occur  as  Wo  decreases,  with 
Kyy  appearing  to  increase  slightly  as  Wo  decreases. 


be  inferred  as  a  residual  term  in  a  heating  budget,  or  it 
must  be  inferred  from  climatological  rainfall  rates.  Un¬ 
fortunately,  latent  heating  is  the  most  important  con¬ 
tribution  to  tropospheric  diabatic  heating  [cf.  Peixoto 
and  Oort ,  1992,  Figure  13.2].  Clearly,  the  number  of 


“reasonable”  modifications  to  the  distribution  shown  in 
Figure  1  is  huge.  We  will  therefore  limit  ourselves  to 
changing  the  strength  of  the  latent  heat  source  Wo  in 
(4).  We  use  a  low  value  of  Wo=2  K  d_1  (experiment  5, 
Table  2)  and  a  high  value  of  Wo=4  Kd"1  (experiment 


Planetary  Wave  Variation 


Figure  8.  (a)  Vertical  wind  averaged  in  time  over  3  model  years  and  horizontally  between  20°S 
and  20° N  <  w*  >20  as  a  function  of  altitude  for  experiments  with  varying  planetary  wave  forcing. 
Long  dashed  curve  shows  <w*f>2o  for  experiment  2  (Znsh= 250  m).  Thin  curve  connecting 

diamonds  shows  <w *  >20  for  experiment  3  {Znsh =125  m).  Dotted  curve  shows  result  for 
experiment  4  {Znsh  =350  m).  (b)  Planetary  wave  EP-flux  divergence  at  30°  N  averaged  over  3 
years  as  a  function  of  altitude.  Line  styles  are  the  same  as  in  Figure  8a.  Additional  dashed  line 
connecting  crosses  shows  Planetary  wave  EP-flux  divergence  for  experiment  1  at  30° N  on  day  16 
of  model  year  3. 
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a) 


a  =  age  jxn-exL+apw  Date  =  3l/12 
b  =  age  jxn-ext2  +  apw  Date=31/12 


c) 


-50  0  50 

Latitude 

c=age  jxn-ext4  +  apw  Date  =  31/12 
d  =  age  jxn-ext  +  apw+ht  Date=3l/12 


Figure  9.  T  as  a  function  of  latitude  and  height  on  day  31  of  model  year  12  for  experiments 
with  fixed  asymmetric  planetary  wave  forcing  amplitudes  and  varying  tropospheric  heating,  (a) 
fto=3  K  d"1  (experiment  1).  (b)  Ho=2  K  d~l  (experiment  5).  (c)  Ho— 4  K  d-1  (experiment  6). 
(d)  Wo=3  K  d""1  with  additional  stratospheric  heat  source  (experiment  7).  Note  that  Figures  4a 
and  9a  are  identical. 


6,  Table  2).  Both  of  these  experiments  use  the  same 
asymmetric  planetary  wave  forcing  and  inputs  to  H  as 
in  Experiment  1  (i.e.,  ZNH=350  m  and  Zsh=175  m). 

Zonal  cross  sections  of  T  at  the  end  of  the  simulations 
using  different  values  of  Ho  (Experiments  1,  5,  and  6) 
are  shown  in  Figures  9a-9c.  A  clear  and  systematic  de¬ 
pendence  on  Ho  is  evident.  There  is  a  nearly  2  year  dif¬ 
ference  in  T  between  the  experiment  with  the  youngest 
middle  atmosphere  and  that  with  the  oldest.  In  ex¬ 
periment  5  with  Ho=l  Kd"1,  the  greater  part  of  the 
middle  atmosphere  above  40  km  is  between  4.5  and  5.0 
years  old  with  portions  of  the  polar  lower  stratosphere 
approaching  6  years  in  age  (Figure  9a).  In  a  20  year 
simulation  with  Wo=2  K  d-1,  the  oldest  air  eventually 
exceeds  6  years  in  age.  This  T  distribution  may  actu¬ 
ally  be  in  better  agreement  with  tracer-based  estimates 
than  our  “baseline”  T  distribution  form  experiment  1. 
In  contrast,  for  experiment  6  with  Wo  =4  K  d_1  (Fig¬ 
ure  9b)  the  bulk  of  the  model  middle  atmosphere  is  less 
than  3.5  years  old.  Planetary  wave  mixing  is  not  re¬ 
sponsible  for  changes  in  T  between  experiments  5  and 
6.  The  Kyy  curves  for  experiments  5  and  6  (Figure  7) 
lie  relatively  close  to  one  another,  with  maximum  NH 
values  between  6xl09  and  8xl09  cm2  s"1.  Given  that 
the  large,  order-of-magnitude  differences  in  Kyy  for  dif¬ 
ferent  Znsh  resulted  in  only  one  year  differences  in  T, 


it  is  unlikely  that  the  small  Kyy  changes  for  different 
Ho  are  producing  the  nearly  2  year  spread  in  middle 
atmospheric  age  between  experiments  5  and  6. 

The  F*  distribution  for  experiments  1,  5,  and  6  are 
shown  in  Figures  10a- 10c.  The  unmixed  age  T*  also 
shows  a  clear  systematic  dependence  on  Ho-  As  Ho  in¬ 
creases  from  2Kd_1to4Kd“1,  the  average  unmixed 
age  at  50  km  decreases  from  around  2.7  years  (Figure 
10b)  to  somewhat  less  than  2.0  years  (Figure  10c).  This 
reduction  in  T*  with  increasing  Ho  is  the  result  of  in¬ 
creased  w*  in  the  tropical  lower  stratosphere,  as  can  be 
seen  in  Figure  11.  Thus  varying  the  imposed  tropical 
tropospheric  heating  has  a  substantial  effect  on  tropi¬ 
cal  ascent  in  our  model,  whereas  varying  the  planetary 
wave  forcing  at  midlatitudes  did  not.  It  is  worth  noting 
that  while  the  unmixed  age  T*  changes  significantly  as 
Ho  varies,  the  changes  in  T*  are  at  most  half  as  large 
as  the  changes  seen  in  the  mixed  age  T.  Thus  a  simple 
increase  in  the  time  required  to  advect  a  parcel  of  air 
into  the  middle  atmosphere  does  not  explain  the  total 
increase  in  T  as  Ho  decreases.  It  appears  that  hori¬ 
zontal  mixing  amplifies  the  delay  in  vertical  transport 
as  the  residual  mean  vertical  velocity  decreases,  even 
though  the  mixing  coefficient  itself  showed  only  a  weak 
dependence  on  Ho- 
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a  =  age*  jxn-ext  +  apw  Date  =  31/l2 
b  =  age*  jxn-ext2+apw  Date  =  3l/12 

Figure  10.  Same  as  Figure  9  except  for 


-50  0  50 

Latitude 


c  =  age*  jxn-ext4  +  apw  Date=31/12 
d  =  age*  jxn-ext+apw+ht  Date=3l/12 

'* .  Note  Figures  4b  and  10a  are  identical. 


4.  Vertical  Motion  in  the  Tropics 

4.1.  Comparison  With  HALOE  Water  Vapor  in 
the  Tropics 

Figure  11  shows  that  the  model’s  vertical  velocity  in 
the  tropical  stratosphere  is  sensitive  to  the  choice  of 
prescribed  tropospheric  heating.  Unfortunately,  direct 
measurements  of  the  residual  mean  vertical  motion  are 
not  available  for  comparison.  However,  recent  satellite 
observations  of  water  vapor  have  given  what  may  be  an 
indirect  observation  of  the  residual  mean  vertical  mo¬ 
tion  in  the  tropical  lower  stratosphere.  The  mixing  ratio 
of  water  vapor  entering  the  stratosphere  in  the  tropics 
is  modulated  by  the  annual  cycle  in  tropical  tropopause 
temperatures  with  cold  tropopause  temperatures  lead¬ 
ing  to  drier  stratospheric  air  [Newell  and  Gould- Stew  art, 
1981;  Mote  et  al ,  1995].  Minimum  tropical  tropopause 
temperatures  during  northern  winter  (<  197  K)  result 
in  water  vapor  mixing  ratios  of  3.5  ppmv.  Water  vapor 
mixing  ratios  rise  to  near  4.5  ppmv  during  northern 
summer  as  tropopause  temperatures  climb  over  201  K. 
This  annually  varying  signal  in  water  vapor  mixing  ratio 
rises  through  the  tropical  lower  stratosphere,  preserving 
its  identity  for  over  a  year  and  a  half  [Mote  et  al.f  1995, 
1996].  This  feature  of  the  tropical  stratosphere  has  been 
recently  dubbed  the  “atmospheric  tape  recorder.” 

Hall  and  Waugh  [1997]  have  used  model  age  spec¬ 
tra  to  argue  that  in  the  tropical  stratosphere  an  an¬ 
nually  varying  signal  like  that  of  water  vapor  will  rise 
with  a  speed  described  by  the  “modal  time”  of  the  age 


spectrum  rather  than  its  first  moment  T.  The  modal 
time  ( tm  in  the  work  of  Hall  and  Waugh  1997)  repre¬ 
sents  the  time  taken  by  the  greatest  fraction  of  mass  to 
reach  a  certain  point,  i.e.,  the  peak  of  the  age  spectrum. 
The  success  of  Lagrangian  parcel  calculations  in  repro¬ 
ducing  the  motion  of  water  vapor  minima  and  maxima 
[e.g.,Mote  et  al .,  1996]  indicates  that  tm  approximates 
the  Lagrangian  travel  time  of  an  air  parcel,  at  least  in 
the  tropical  lower  stratosphere.  In  other  words,  in  the 
tropics  we  should  have  tm& r*,  and  we  therefore  ex¬ 
pect  the  vertical  motion  of  the  unmixed  age  tracer  /i* 
to  compare  well  with  tropical  water  vapor  observations. 

Plates  la-lc  show  the  evolution  of  latitude- averaged 
unmixed  age  tracer  </i*  >12  as  a  function  of  height  and 
time  for  the  experiments  with  varying  Wo  (experiments 
1,5,  and  6).  Below  18  km  the  tracers  are  rapidly  mixed 
and/or  advected  vertically  so  that  little  or  no  vertical 
gradient  is  evident  in  the  isopleths  of  <  /i*  >12.  Be¬ 
tween  18  km  and  about  35  km,  a  region  of  slow  upward 
drift  exists.  In  this  layer  the  sensitivity  of  vertical  mo¬ 
tion  to  Wo  is  clearly  evident.  By  35  km  the  isopleths  of 
<  /i*  >12  from  Experiment  6  with  Wo  =4  K  d_1  (Plate 
lc)  are  about  0.75  years  ahead  of  those  from  Experi¬ 
ment  5  with  Wo=2  K  d_1  (Plate  la).  Above  35  km 
the  isopleths  from  the  different  experiments  again  be¬ 
come  steeper  and  nearly  parallel,  reflecting  the  increas¬ 
ingly  vigorous  meridional  circulation  at  these  altitudes. 
Thus,  most  of  the  difference  in  unmixed  age  of  air  be¬ 
tween  experiments  develops  in  the  slow  ascent  region 
between  18  and  35  km.  The  same  is  true  for  mixed 
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Figure  11.  Vertical  wind  averaged  in  time  over  sev¬ 
eral  model  years  and  horizontally  between  20° S  and 
20°N  (Cw*’*  >20)  as  a  function  of  altitude  for  exper¬ 
iments  with  varying  tropospheric  heating.  Solid  curve 
shows  <w*t  >20  for  experiment  1  (Ho =3  Kd*1).  Long- 
dashed  curve  shows  <tu**>2o  for  experiment  5  (Ho=2 
K  d-1).  Dotted  curve  shows  <w*f> 20  for  experiment 
6  (Hq=4  K  d"1).  Thin-dashed  lines  connecting  open 
diamonds  show  <w**>20  for  experiment  7  (Hq=3  K 

d'1  with  Tiatrat )•  Note  that  <^**>20  for  experiments 
1  and  7  converge  below  15  km.  Dashed  line  joining  solid 

circles  shows  a  profile  of  <H/QZ  >20  j  an  approxima¬ 
tion  to  the  steady  state  vertical  velocity,  for  experiment 
1. 


age,  which  shows  a  qualititavely  similar  time  evolution 
in  the  tropics,  except  that  the  isopleths  of  mixed  age 
rise  more  slowly. 

The  color  background  for  Plates  la-lc  shows  HALOE 
2(CH4)+H20  mixing  ratio  data  from  the  tropics  aver¬ 
aged  between  8.5°S  and  8.5°N  as  a  function  of  height 
and  time.  Similar  figures  based  on  MLS  as  well  as 
HALOE  H2O  measurements  have  been  shown  elsewhere 
[e.g.,  Mote  et  a/.,  1995,  1996].  We  show  them  in  Plate  1 
simply  to  facilitate  comparison  with  the  present  model 
results.  The  overall  pattern  of  time  evolution  of  the 
model  <//*  >12  is  very  similar  to  that  of  the  HALOE 
data.  The  HALOE  data  in  Plate  1  exhibit  a  layer  of 
slow  ascent  between  18  and  30  km.  This  is  evident  in 
the  shallow  slope  of  minima  and  maxima  of  water  vapor 


or  2(CH4)+H2  0  in  time-height  plane.  Above  30  km  the 
“tape  recorder”  signal  weakens  quickly,  but  it  appears 
that  ascent  is  becoming  more  rapid.  Below  20  km  the 
data  also  suggest  more  rapid  vertical  transport  (this  is 
clearer  in  Plates  la  and  lb  of  Mote  et  al.  [1996]).  The 
traces  of  the  water  vapor  minima  and  maxima  in  the 
HALOE  data  imply  a  typical  transit  time  of  about  1.0 
year  between  20  km  and  30  km,  although  there  is  signif¬ 
icant  interannual  variability.  Note  the  relatively  rapid 
ascent  between  July  1993  and  April  1994  and  the  rel¬ 
atively  slow  ascent  commencing  in  July  1995.  Overall, 
there  is  good  quantitative  agreement  between  <//*  >12 
in  experiments  1  and  5  in  the  18  to  35  km  layer  and 
HALOE  2(CH4)+H20,  while  ascent  in  experiment  6 
(Ho=4  K  d-1)  is  somewhat  too  rapid. 

4.2.  Forcing  of  Model  Vertical  Motion  in  the 
Tropics 

The  sensitivity  of  w*  and  T  to  prescribed  tropospheric 
heating  shows  that  poorly  or  incompletely  simulated 
tropospheric  processes  in  a  2-D  zonally  averaged  model 
can  have  an  impact  well  into  the  middle  atmosphere. 
Identifying  the  ultimate  forcing  for  vertical  motion  in 
the  tropical  stratosphere  in  a  nonlinear  2-D  model  is 
complicated  by  the  fact  that  both  drag  and  heating  pro¬ 
duce  nonlocal  circulations,  which  then  lead  to  nonlocal 
changes  in  temperature  and  zonal  wind.  Changes  in 
zonal  wind  or  temperature  can  of  course  lead  to  further 
changes  in  drag  and  heating  and  so  on.  Remote  forc¬ 
ing  by  midlatitude  drag  is  currently  believed  to  be  the 
most  important  forcing  for  vertical  motion  in  the  trop¬ 
ical  lower  stratosphere  [e.g.,  Iwasaki ,  1992;  Yulaeva  et 
al,  1994;  Holton  et  al.,  1995].  However,  in  this  section 
we  will  show  that  the  sensitivity  of  our  age  simulations 
to  tropical  heating  occurs  without  a  significant  contri¬ 
bution  from  midlatitude  planetary  wave  drag.  First, 
we  examine  in  more  detail  the  time  behavior  of  w*  and 
components  (see  below)  of  w*  in  the  tropical  strato¬ 
sphere  during  the  age  simulations  in  section  3.  Then 
we  will  examine  the  changes  to  the  stream  function  as 
Ho  varies. 

4.2.1.  Decomposition  of  the  residual  circula¬ 
tion.  .  The  residual  circulation  in  the  quasigeostrophic 
TEM  formulation  is  obtained  from  a  linear  elliptic  equa¬ 
tion  forced  by  gradients  of  the  momentum  forcing  X 
(2)  and  total  heating  H  (3),  [e.g.,  Andrews  et  al  1987, 
Equation  3.5.7].  The  total  residual  circulation  can  thus 
be  broken  down  into  a  sum  of  separate  circulations 
forced  by  the  individual  terms,  or  groups  of  terms,  in 
X  and  H.  Each  of  these  component  circulations  is 
that  which  would  result  if  no  additional  atmospheric 
response  occurred  as  a  result  of  the  particular  compo¬ 
nent  circulation  under  consideration  (e.g.,  the  adiabatic 
responses  to  momentum  forcing  shown  by  Haynes  et  al. 
[1991]  and  by  Holton  et  al.  [1995,  Figure  4a].  The 
Appendix  analyzes  the  development  of  “downward  con¬ 
trol”  [Haynes  et  a/.,  1991]  in  terms  of  such  component 
stream  functions. 

4.2.2.  Seasonal  variation  of  vertical  velocity 
in  the  tropics  .  Figure  12  shows  model  time  series 


UNMIXED  Age  Tracer  jx 
superimposed  on  HALOE  2*CH4+H20 


11,278 


BACMEISTER  ET  AL.:  AGE  OF  AIR  IN  A  ZONALLY  AVERAGED  MODEL 


"5  v'O  g  <b  .a  -A 

.a  S-a  I  %  J  | 

£  ~  * 

C  73  CO  <U 

p  fl  £  <D  o 

8  *5  m  is 


o  *.  ^ 

_  ^  n  ®  h  ■*»  - 

•2  .SP^  |  toe- 

2  -S9.  j 
:s  gtjgj-g'i? 

a--g5§«  si 

bp  fl  W  +3  *g  «  W 

a  w  s?  ^  , 

£  <e2  g.«  5 

ir  v  i— I  bO-u  ,-g 

*®  §  ^  -3  §  "a  § 

T3  .^SEg 

-a  a  8  ^  s  ^  s 
o*  -ag  g\s  g< 
^  ^  P  05  w  s  v 


y  T-H  ^  5 

S  O 


« a  a 

■'*  g.  g 

II  O.’B 

o  td  fl 
>*  O 

sS-S  ° 

CO  13  XJ 

'T7  fll 

a 

III 
11 S 

■  X  *  -d 
$ 

£ 

X§  o 
7  «  >> 

*1  § 

CN  § 

II  8  ^ 
o^d  '■d 
* 


»  <U 
|-H  HS  ' 

.  o 

°  6  - 
co  G 


>£  * 


i-o  ^ 

-O  Pi 
fl  P 


ouco 

■|H 


£  8 

AS  CN 


*  £ 
>h  <D 

s  >» 

a; 

aw 

gO 

*  W 

3« 


^  fl  *fl  ^ 

« g>|s  s 

,fh  .5  ^  ^  ^ 
2.3--S® 

o  3  ^  y 

£  SP 
a  2  ° 

*h  aj 

a  a) 


>>  a> 

®  -4f 

SO?  fl 
^  f1  CJ 
O  %  rj 

^s-|gas|“ 

o  ^  ^2  fl  Sfl  o  hn 

civ  fl  .fl  Mi  *h  +5  0)  wO 

0)  5P  fl  ^  Oh  <U  So^fl 

|*§  2  go  5 ST- fa 
ce>ee^5K..oflo 

•  r  A  *■  fl  A  o 

fl  "  ' 


fl  >  —  . — . 
b  S  a  o 
O  u  g  O 

fJL  fl  O  tfc  >-— 1 ' — -  >-r  i— 

8  s-s^-gl 

^  -fix*  ^  §  7!  <R  5 


o  ^  *  fl 

-  60  Oh 


b  -fl  -r  .  fl  co  ^ 

O  boO  2  <L)  •  O  o 
w  -Si  «  <d  o  lo  cm  rr 

I 

S?  5 


g  fl  IN  «  bo-o  |  W  -  & 
•2  -  H  g-.B  5  2  S  « 

S  O  ^fl  ^  .  A  ^ 

>  <u 3  §2^  ^ ^  bo^r s  ^ 
r  «>  «  g  (g  «  flEc# 

Mh  fl  O  C n  tw  -P  fl  fl  H  A4  oq  boW  +->  bO'» _ ^o 


7Z»  tn  -*-3  +3 
«  <L>  <U  Oh  ^  fl 

T  up  >  «  o  3  « 


LO 


O  fl 

OQ  o 


fl  ^  ^ 
fl  -4-3  ^  fl 

•  A.y  _ 
2f£  £ 


co  n 
a)  fl; 
A  Oh 


<D  ’ 


BACMEISTER  ET  AL.:  AGE  OF  AIR  IN  A  ZONALLY  AVERAGED  MODEL 


11,279 


Figure  12.  Time  series  of  various  components  of  tropically  averaged  residual  vertical  velocity 
<w*  >20  at  18  km  (bottom)  and  29  km  (top)  for  experiment  1.  The  solid  curves  show  the  total 
vertical  velocity  <  w*  >20-  Short-dashed  curves  show  the  drag  induced  component  <  w %  >20,  and 
long-dashed  curves  show  the  thermally  induced  component  <w ^  >20 •  Circles  show  a  time-series 
of  <  Ti/Qz  >20  j  an  approximation  to  the  steady  state  vertical  velocity. 


of  various  components  of  <  w*  >20  for  experiment  1  at 
18  km  and  29  km.  The  vertical  velocity  has  been  bro¬ 
ken  down  into  two  components:  which  is  forced  by 

the  total  drag,  and  which  is  forced  by  the  total  net 
heating.  The  sum  of  these  is  the  total  residual  vertical 
velocity  w* .  Figure  12  shows  that  <  w*  >20  as  well  as 
<  u>x  >20  and  <  w ^  >20  have  a  clear  annual  cycle  at  18 
km.  The  amplitude  of  the  annual  cycle  in  w*  is  around 
0.005  cm  s"1,  which  yields  a  roughly  estimated  heat¬ 
ing/cooling  rate  of  0.05  K  d-1,  more  than  sufficient  to 
account  for  the  model’s  1-2  K  annual  cycle  at  the  trop¬ 
ical  tropopause.  However,  in  our  model  this  annually 
varying  signal  is  superimposed  on  a  “floor”  value  of  w* 
close  to  0.020  cm  s'1  as  indicated  in  the  figure.  Thus 
in  our  model  the  amplitude  of  the  annually  varying  sig¬ 
nal  in  w*  is  around  one  third  of  this  floor  value.  This 
means  that  although  significant  wave-induced  annual 
variations  in  tropopause  w*  occur  in  our  model,  they 
are  superimposed  on  a  relatively  large  time  indepen¬ 
dent  flux  through  the  tropopause.  The  total  amplitude 
of  the  annual  cycle  of  tropical  w*  at  18  km  comes  in 
roughly  equal  parts  from  w %  and  w However,  the 
“DC”  portion  of  the  time  series  of  <  w*  >20  is  almost  en¬ 
tirely  contributed  by  w The  phasing  of  <  w %  >20  and 
<^>20  is  slightly  different,  with  w'x  peaking  around 
January  1  and  w ^  peaking  in  March.  The  presence  of 


an  annual  cycle  in  w ^  is  a  clear  indication  that  midlat¬ 
itude  drag  is  indirectly  inducing  a  thermally  forced  cir¬ 
culation,  as  discussed_in  the  Appendix.  Note  also  that 

<  w h  >20  and  <  Ti/Qz  >20  are  slightly  out  of  phase, 
which  is  another  indication  of  remote  forcing. 

At  29  km  the  situation  is  quite  different.  Little  an¬ 
nual  cycle  is  evident,  but  a  strong  semi-annual  cycle 
is  visible  in  all  components  of  <  w*  >20*  The  mo¬ 
mentum  forced  circulation  w'x  contributes  a  smaller 
fraction  to  the  semiannual  cycle  in  tropical  w*  than 
does  Wft*  The  phase  of  <  w^  >20  at  29  km  is  close 
to  that  of  <  /H/QZ  >20  with  maxima  occurring  near 
both  equinoxes.  In  contrast,  the  phase  of  <  w'x  >20  is 
shifted  such  that  maxima  occur  somewhat  close  to  the 
solstices.  As  was  the  case  at  18  km,  the  mean  value  of 

<  w*  >20  is  substantially  larger  than  the  amplitude  of 
the  time-varying  portion  and  is  contributed  primarily 
by  <  w^  >20-  The  weak  annual  cycle  and  strong  semian¬ 
nual  cycle  in  <w*  >20  at  29  km  is  consistent  with  trop¬ 
ical  or  hemispherically  symmetric  forcing  rather  than 
with  asymmetric  midlatitude  planetary  wave  drag.  The 
relationship  of  our  models  SAO  in  zonal  wind,  which  is 
weak  below  30  km  [Summers  et  al.}  1997],  to  the  SAO 
in  w*  is  not  yet  clear. 

The  presence  of  annual  cycle  in  <  w*  >20  at  18  km 
agrees  with  the  idea  that  vertical  motion  through  the 
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Figure  13.  Three  year  average  of  mass  stream  function  pip1  and  components  for  experiment 
1  (W0  =3  K  d-1):  (a)  pip*  total  mass  streamfunction;  (b)  pipx  mass  streamfunction  forced  by 
total  mechanical  forcing;  (c)  pipn  mass  streamfunction  forced  by  total  thermodynamic  forcing; 
(d)  pippw*  mass  streamfunction  forced  by  drag  from  planetary  waves.  Solid  contours  are  positive 
and  indicate  a  clockwise  circulation.  Thick  contours  are  plotted  every  2.5  xlO7  mbar  cm2  s”1, 
while  thin  contours  are  plotted  every  2.5  xlO6  mbar  cm2  s-1. 


tropical  tropopause  is  modulated  by  extratropical  drag 
[e.g.,  Yulaeva  et  al. ,  1994].  However,  even  at  18  km  in 
our  model  the  annually  varying  component  of  w*  repre¬ 
sents  at  most  one  fifth  of  the  time-averaged  w*  at  the 
tropical  tropopause.  It  should  be  pointed  out  again  that 
observations  indicate  a  larger  annual  cycle  in  tropical 
tropopause  temperature  than  is  present  in  our  model,  so 
in  reality,  a  somewhat  larger  fraction  of  the  annually  av¬ 
eraged  transport  through  the  tropical  tropopause  may 
occur  as  the  result  of  extratropical  driving.  In  contrast, 
at  29  km  both  observations  and  model  results  indicate 
a  weak  annual  cycle  (<1  K)  in  temperature.  Thus  there 
is  no  a  priori  reason  to  invoke  extratropical  control  of 
vertical  transport  at  these  altitudes. 

4.2.3.  Time-averaged  stream  function.  From 
Figure  12  it  is  clear  that  in  our  model  the  dominant 
mode  of  vertical  transport  throughout  the  tropical  strato¬ 
sphere  is  advection  by  the  time-mean  portion  of  the 
vertical  velocity.  Both  annual  and  semiannual  varia¬ 
tions  in  the  model  are  small  compared  to  time-mean 
value  <w*t>20-  Figure  13  shows  plots  of  multiyear 


averages  of  the  total  meridional  mass  stream  function 
pip*  (Figure  13a)  as  well  as  several  of  its  components 
for  Experiment  1.  The  mass  stream  function  forced  by 
total  drag  p ip^  (Figure  13b)  is  shown  along  with  that 

forced  by  total  net  heating  pip^  (Figure  13c).  Consis¬ 
tent  with  the  time  series  of  w* ,  ,  and  w x  in  Figure 

12  the  time-mean,  total  mass  stream  function  pip*  in 
the  tropics  is  dominated  by  pip^  .  Also  shown  is  the 
subcomponent  of  pipx*  driven  by  only  planetary  wave 
EP-flux  divergence  -pV'pw*  (Figure  13d).  Strong  hemi¬ 
spheric  asymmetry  is  evident  in  both  pip*x  and  pt/>pW 
as  is  expected  for  an  experiment  with  asymmetric  plan¬ 
etary  wave  forcing.  However,  pip^  shows  little  hemi¬ 
spheric  asymmetry,  and  above  20  km  is  more  concen¬ 
trated  in  the  tropics  than  the  drag-induced  circulations. 
This  argues  against  a  large  remote  effect  by  extratrop¬ 
ical  planetary  wave  drag  on  the  time-mean  net  heating 
of  the  tropical  stratosphere  in  the  present  model. 

Figure  14  shows  the  difference  in  these  mass  stream 
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Figure  14.  As  in  Figure  13  except  that  the  difference  in  3  year  average  mass  streamfunc- 
tion  components  between  experiments  1  (7f0=3  K  d"1)  and  6  ( Hq-A  K  d"1)  is  shown.  We 

show  averaged  components  of  A  pip*  =  pip*  Exp  6  —  pip*  Exp  j .  so  positive  values  indicate  stronger 
circulations  in  experiment  6. 


function  components  between  experiment  1  and  Exper¬ 
iment  6.  As  the  prescribed  tropospheric  heating  is  in¬ 
creased  the  change  in  mass  stream  function  follows  the 
pattern  in  Figure  13;  the  momentum  forced  components 
change  asymmetrically,  while  the  thermally  forced  com¬ 
ponent  changes  symmetrically  about  the  Equator.  Fur¬ 
thermore,  the  changes  in  (Figure  14b)  and 
(Figure  14d),  while  significant,  are  confined  in  the  verti¬ 
cal  and  thus  cannot  explain  the  changes  in  w**  shown  in 
Figure  11,  which  extend  over  a  deep  layer  of  the  tropi¬ 
cal  stratosphere.  It  is  these  changes  in  w*  which  lead  to 
the  sensivity  of  T  to  the  prescribed  tropospheric  heat¬ 
ing.  Similar  changes  to  p and  its  components  are 
observed  when  going  to  reduced  tropospheric  heating, 
i.e.,  between  experiment  1  and  experiment  5. 

From  the  results  shown  in  Figures  12-14  we  con¬ 
clude  that  tropical  ascent  in  our  model  is  controlled 
by  tropical  or  at  least  hemispherically  symmetric  pro¬ 
cesses.  Near  tropopause  altitudes  the  model  may  under¬ 
estimate  the  influence  of  midlatitude  planetary  waves. 
However,  if  we  assume  the  relative  importance  of  hemi¬ 
spherically  asymmetric  planetary  wave  drag  at  the  trop¬ 
ical  tropopause  can  be  estimated  by  comparing  the  am¬ 


plitude  of  the  annual  cycle  in  w*  to  the  background 
floor  value  of  w*  at  18  km  in  Figure  12,  then  even  dou¬ 
bling  the  planetary  wave  effect  will  only  result  in  an 
annually  varying  component  that  is  at  most  of  similar 
amplitude  to  the  floor  value.  Doubling  the  amplitude 
of  the  annual  oscillation  in  w*  at  18  km  would  also  dou¬ 
ble  the  amplitude  of  the  model’s  annual  oscillation  in 
tropical  tropopause  temperature,  which  is  small  (3  K) 
but  within  the  range  of  observed  amplitudes  (Section 
3.1).  Doubling  the  amplitude  of  the  model’s  annual  cy¬ 
cle  in  tropical  tropopause  temperature  would  place  it 
well  within  if  not  at  the  high  end  of  the  observations. 
Thus  we  conclude  further  that  although  our  model  may 
underestimate  the  contribution  of  planetary  wave  drag 
on  tropical  ascent  at  the  tropopause,  it  does  so  by  less 
than  a  factor  of  2.  At  higher  altitudes  in  tropical  strato¬ 
sphere  we  find  no  compelling  reason  to  disbelieve  the 
model’s  simulation  of  vertical  motion. 

4.3.  Response  to  Lower  Stratospheric  Heat 
Source 

Some  of  the  difficulty  in  identifying  the  cause  of  the 
model’s  sensitivity  to  variations  in  tropospheric  heating 
arises  from  the  fact  we  expect  a  strong  response  in  the 
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Figure  15.  (a)  Differences  in  net  heating  averaged  over  3  model  years  and  between  20°S 

and  20° N  <  Tt  >20  as  a  function  of  altitude.  Solid  curve  shows  net  heating  change  between 
experiments  6  and  1,  <  >20  Exp  6  -  <  W_> 20  Exp  i-  Dashed  curve  shows  change  between 
experiments  7  and  1,  <  W  >20  Exp  7  —  <  W  >20  Exp  l-  (b)  RMS  changes  to  temperature 
averaged  over  3  model  years  and  between  20°S  and  20°N.  Solid  curve  shows  <  AT2  >\^  for 
experiments  6  and  1,  dashed  curve  shows  <  AT2  >\^  for  experiments  7  and  1.  In  both  instances 
the  large  values  of  <  AT2  >2q2  correspond  to  temperature  increases  in  the  vicinity  of  additional 
imposed  heating,  i.e.,  more  imposed  tropospheric  heating  in  experiment  6  and  more  imposed 
stratospheric  heating  in  experiment  7. 


tropospheric  Hadley  circulation  to  such  variations.  The 
Hadley  cell  response  can  lead  to  changes  in  low-level 
zonal  winds  as  well  as  in  planetary  wave  and  gravity 
wave  generation,  which  can  in  turn  modify  the  momen¬ 
tum  forcing  of  the  meridional  circulation.  The  existence 
of  such  indirect  effects  on  the  model’s  atmospheric  drag 
field  is  easily  seen  in  Figure  14d  where  changes  in  pre¬ 
scribed  tropospheric  heating  led  to  differences  in  the  cir¬ 
culation  forced  by  planetary  wave  EP-flux  divergence. 
The  changes  in  <^**>20  with  increasing  Wo  in  Fig¬ 
ure  11  are  associated  with  small  changes  in  the  strato¬ 
spheric  net  heating.  It  is  of  interest  to  see  whether 
the  circulation  in  the  model  is  sensitive  to  direct  addi¬ 
tion  of  small  amounts  of  heat  to  the  lower  stratosphere. 
This  should  leave  the  model’s  Hadley  circulation  rela¬ 
tively  unaffected,  while  still  producing  changes  to  the 
stratospheric  net  heating  similar  to  those  resulting  from 
changes  to  Wo-  Furthermore,  such  an  experiment  ad¬ 
dresses  real  uncertainties  in  stratospheric  heating  which 
result,  for  example,  from  volcanic  aerosols  [e.g.,  Kinne 
et  ai,  1992;  Kinnison  et  al,  1994;  Rosenfield  et  al, 
1997]. 

We  perform  a  final  experiment  (experiment  7)  us¬ 
ing  the  same  value  of  Wo  as  in  experiment  1,  but  with 
a  weak  additional  heat  source  centered  in  the  tropical 
lower  stratosphere  as  shown  in  Figure  lb.  The  magni¬ 
tude  of  this  heat  source,  0.3  K  d”1,  was  chosen  so  that 
the  final  net  heating  difference  in  the  lower  stratosphere 
was  similar  to  that  obtained  between  experiment  1  and 


experiment  6  by  increasing  Wo-  Figure  15a  shows  the 
difference  in  tropical  net  heating  profiles  <  W  >20  be¬ 
tween  experiments  1  and  6  (solid  line)  and  between  ex¬ 
periments  1  and  7  (dashed  line).  It  is  clear  from  Figure 
15  that  net  heating  changes  of  0.05  to  0.1  K  d_1  are  ob¬ 
tained  in  a  deep  layer  of  the  tropical  stratosphere  both 
by  increasing  tropospheric  heating,  as  in  going  from  ex¬ 
periment  1  to  experiment  6,  or  by  adding  heat  directly 
to  the  lower  stratosphere,  as  in  experiment  7.  Inter¬ 
estingly,  the  stratospheric  temperature  change  between 
experiments  1  and  6  (Figure  15b)  is  small,  implying  that 
any  increase  in  thermal  forcing  going  from  experiments 
1  to  6  is  nearly  balanced  by  increased  w* ,  at  least  in 
the  stratosphere.  In  contrast,  the  tropical  stratosphere 
is  2-3  K  warmer  in  experiment  7  than  in  experiment 
1.  This  indicates  that  both  increased  IR  cooling  and 
increased  w*  occur  when  heat  is  added  directly  to  the 
stratosphere. 

Despite  the  differences  in  stratospheric  temperature 
profiles,  the  stratospheric  circulations  in  experiments  6 
and  7  are  quite  similar.  Figures  9d  and  lOd  show  T  and 
r*  for  experiment  7  after  11  model  years.  There  are 
clear  similarities  between  these  and  the  T  and  T*  dis¬ 
tributions  for  experiment  6  (Figures  9c  and  10c),  con¬ 
firming  that  the  restricted  stratospheric  heat  source  has 
an  overall  effect  on  the  circulation  equivalent  to  that  of 
simply  raising  Wo  • 

As  expected,  the  addition  of  a  tropical  stratospheric 
heat  source  did  not  affect  midlatitude  drag  to  the  same 
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Figure  16.  As  in  Figure  14  except  for  experiments  1  (%o=3  K  d”1)  and  7  (7io=3  K  d_1  with 
Wstrat).  We  show  averaged  components  of  Apip*f  =  pip*^^  6  -  pip* gx 


degree  as  an  increase  in  Wo>  This  can  be  seen  in  Fig¬ 
ure  16  which  shows  the  change  in  components  of  the 
mass  stream  function  pip*  between  experiments  1  and 
7.  The  stratospheric  changes  to  the  total  stream  func¬ 
tion  pip*  and  to  the  thermally  forced  component  pip ^ 
(Figures  16a  and  16c)  are  of  similar  magnitude  and 
shape  to  those  obtained  between  Experiments  1  and 
6  (Figures  14a  and  14c).  Note  that  the  changes  in  pip* 
and  pipft  in  Figures  16a  and  16c  die  off  rapidly  be¬ 
low  15  km.  The  drag  induced  circulations  pip *x%  and 

p^PW  show  much  weaker  changes  overall  between  Ex¬ 
periments  1  and  7  (Figures  16b  and  16d)  than  were 
obtained  between  Experiments  1  and  6  (Figures  14b 
and  14d).  Thus,  the  addition  of  a  stratospheric  heat 
source  intensifies  tropical  vertical  motion  (and  reduces 
T)  without  creating  a  significant  change  to  the  circula¬ 
tion  forced  by  midlatitude  drag. 

The  sensitivity  of  our  model's  circulation  and  T  to 
tropospheric  heating  variations  must  be  viewed  with 
caution  due  to  the  incomplete  tropospheric  dynamics 
in  the  model.  However,  the  response  of  the  model  to 
the  restricted  stratospheric  heat  source  given  by  (4d) 
was  obtained  without  appreciable  effects  in  the  model 
troposphere.  The  ~0. 10  K  d_1  net  heating  change  in¬ 


duced  by  the  stratospheric  heat  source  (Figure  15a)  is 
also  well  within  the  possible  perturbations  induced  by 
a  major  volcanic  eruption. 

We  have  not  yet  isolated  the  mechanism  responsible 
for  increasing  our  model’s  tropical  vertical  velocities  as 
Wo  increases.  However,  the  analysis  shown  here  indi¬ 
cates  that  (1)  midlatitude  wave  drag  is  not  involved, 
and  (2)  similar  effects  on  w*  can  be  induced  by  a  purely 
stratospheric  heat  source.  There  are  at  least  two  pos¬ 
sibilities  which  have  not  been  addressed  here.  (1)  Tro¬ 
pospheric  temperature  increases  as  Wo  grows  (Figure 
15b),  are  “seen”  by  the  stratosphere  and  produce  lo¬ 
cal  increases  in  stratospheric  net  heating.  (2)  Increased 
w *  in  the  lower  stratosphere  is  remotely  forced  by  in¬ 
creased  Wiatent-  Either  of  these  mechanisms  could  be 
reinforced  by  secondary  circulations  arising  from  the 
initial  changes  to  the  lower  stratospheric  net  heating. 
A  complete  examination  of  these  possibilities  is  beyond 
the  scope  of  the  present  paper. 

5.  Summary  and  Discussion 

This  study  examined  factors  controlling  age  of  air  in 
a  2-D,  zonally  averaged  model  of  the  atmosphere.  We 
examined  age  distributions  derived  from  two  tracers. 
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The  forcing  and  evolution  of  both  tracers  were  identi¬ 
cal  in  each  model  simulation  except  that  one  tracer  was 
affected  by  horizontal  mixing  while  the  other  was  not. 
Comparisons  of  mixed  age  (r)  and  unmixed  age  (P) 
in  our  simulations  reveal  qualitative  as  well  as  quanti¬ 
tative  differences.  In  all  of  the  simulations,  P  was  1  to 
2  years  less  than  T  in  most  of  the  middle  atmosphere. 
The  polar  lower  stratosphere  was  the  only  region  of  the 
atmosphere  in  which  V*  was  larger  than  T.  We  found 
a  systematic  sensitivity  of  P  to  the  strength  of  pre¬ 
scribed  tropical  tropospheric  heating  in  the  model.  As 
the  peak  strength  of  the  prescribed  heating  increased 
from  2  K  d"1  to  4  K  d"1,  with  fixed  planetary  wave 
forcing,  the  typical  middle  atmospheric  value  of  P  de¬ 
creased  from  around  2.5  to  less  than  1.7  years  (Figure 
10).  In  contrast,  when  prescribed  tropospheric  heating 
was  fixed  and  planetary  wave  forcing  amplitudes  were 
varied  from  125  m  to  350  m,  little  change  was  noted  in 
P  (Figure  6).  Thus,  T*  (and  w*)  in  our  model  was  sen¬ 
sitive  to  changes  in  prescribed  tropospheric  heating  but 
not  to  changes  in  planetary  wave  forcing  amplitudes. 

Mixed  age  T  was  found  to  be  sensitive  to  changes, 
in  both  prescribed  tropospheric  heating  and  planetary 
wave  forcing  amplitudes.  As  planetary  wave  forcing  am¬ 
plitudes  were  increased  from  125  to  350  m,  Kyy  in  the 
lower  stratosphere  increased  by  a  factor  of  8  (Figure  8). 
This  stronger  mixing  resulted  in  increased  values  of  T  in 
the  middle  atmosphere  for  larger  planetary  wave  forc¬ 
ing  even  though  P  stayed  relatively  constant  (Figure 
6).  Stronger  vertical  motion  in  the  tropics  also  led  to 
lower  values  of  T  as  Wo  increased  (Figure  9).  A  compar¬ 
ison  of  T  with  atmospheric  observations  yielded  reason¬ 
able  agreement  for  Wo  =2  or  3  K  d"-1.  Our  peak  lower 
stratospheric  age  of  ~6  years  in  experiment  5  (Wo =2  K 
d-1)  agrees  with  age  estimates  based  on  balloon-borne 
CO2  [Schmidt  and  Khedim ,  1991]  but  is  substantially 
less  than  the  very  high  values  of  8  to  10  years  obtained 
by  Harnisch  et  ai  [1996]  from  their  balloonborne  SF6 
measurements.  Our  upper  stratospheric  T  values  of  4 
to  5  years  for  Wo =2  and  3  K  d"1  were  also  somewhat 
below  the  5.9  year  age  estimate  from  HALOE  HF*  [Rus- 
sell  et  a/.,  1996].  Finally,  using  Wo=2  Kd"1,  we  obtain 
reasonable  agreement  with  the  3-D  MACCM2  age  simu¬ 
lations  of  Waugh  et  ai  [1997]  using  the  updated  gravity 
wave  parameterization  (MACCM2  version  2). 

The  high  values  of  midst r at ospheric  T  (>5  years)  ob¬ 
tained  from  gases  with  linearly  increasing  tropospheric 
sources  (e.g.,  HF,  C02,  SF6)  are  much  longer  than  the 
transit  times  of  tropical  water  vapor  pulses  in  HALOE 
and  MLS  observations  of  (i.e.,  the  atmospheric  “tape 
recorder”).  These  pulses  travel  from  15  km  to  35  km  in 
approximately  1.2  -  1.5  years  (see  Plate  1;  or  Mote  et 
ai,  [1996,  Plate  1]).  Mote  et  ai  [1996]  show  that  the 
minima  and  maxima  of  zonally  averaged  tropical  wa¬ 
ter  vapor  approximately  trace  the  path  described  by  a 
Lagrangian  parcel  in  the  mean  meridional  circulation. 
The  transit  time  of  such  pulses  may  be  a  better  ap¬ 
proximation  of  the  modal  time  (or  peak)  of  the  age 
spectrum  than  its  first  moment  T,  which  corresponds  to 
“age”  derived  from  long-lived  tracers  [Hall  and  Waugh , 


1997].  The  time  evolution  of  the  unmixed  age  tracer 
fi*  in  the  tropics  (Plate  1)  shows  that  with  Wo=2  K 
d-1  and  Wo=3  K  d-1,  our  model  does  a  good  job  of 
reproducing  the  apparent  vertical  trajectories  of  pulses 
in  2CH4+H2O  observed  by  HALOE.  Thus  simulations 
with  these  values  of  Wo  quantitatively  reproduce  both 
the  high  values  of  T  estimated  for  the  midstratosphere 
using  trace  gas  mixing  ratios  and  the  apparent  config¬ 
uration  of  tropical  ascent  inferred  from  HALOE  and 
MLS  water  vapor  observations.  This  indicates  that  our 
model  correctly  simulates  the  balance  between  verti¬ 
cal  advection  and  horizontal  tropical/midlatitude  ex¬ 
change  which  characterizes  tracer  transport  in  the  trop¬ 
ical  stratosphere  [e.g.,  Plumb)  1996;  Volk  et  a/.,  1996; 
Minschwaner  et  a/.,  1996].  Put  another  way,  these  re¬ 
sults  show  that  our  simulations  with  Wo=2  K  d-1  and 
Wo =3  K  d-1  reproduce  both  the  first  moment  and  the 
peak  of  the  atmospheric  age  spectrum.  It  is  gratifying 
that  a  correct  simulation  of  these  processes  is  associated 
with  planetary  wave  and  tropospheric  heating  param¬ 
eters  that  also  yield  reasonable  simulations  of  strato¬ 
spheric  zonal  winds  and  temperatures  in  the  model. 

The  sensitivity  of  T  and  T*  in  our  model  to  prescribed 
tropical  tropospheric  heating  Wo  is  due  to  vertical  ve¬ 
locities  in  the  tropical  stratosphere,  which  exhibited  a 
larger  than  expected  senstivity  to  Wo,  while  exhibit¬ 
ing  little  if  any  sensitivity  to  changes  in  planetary  wave 
forcing  amplitudes.  Changes  to  Wo  can  affect  the  global 
distribution  of  either  drag  or  net  heating,  both  of  which 
could  impact  the  meridional  circulation.  By  decom¬ 
posing  the  total  residual  stream  function  into  separate 
components  associated  with  drag  X  and  net  heating  W, 
we  showed  that  changes  in  the  tropical  vertical  velocity 
are  not  forced  by  changes  to  the  model  planetary  wave 
drag.  We  found  that  the  total  residual  stream  function 
above  25  km  in  the  tropics  was  dominated  by  the  con¬ 
tribution  from  the  heating.  Furthermore,  even  though 
hemispherically  asymmetric  planetary  wave  forcing  was 
used,  we  found  that  the  change  to  the  residual  stream 
function  as  Wo  varied  remained  relatively  symmetric. 
These  facts  indicate  that  varying  the  prescribed  tropo¬ 
spheric  heating  in  our  model  alters  tropical  ascent  pri¬ 
marily  by  an  upward  directed  change  to  the  net  heat¬ 
ing  of  the  tropical  lower  stratosphere,  rather  than  by  an 
indirect,  mechanically-forced  path  mediated  by  midlat¬ 
itude  drag. 

The  tropospheric  response  of  our  2-D,  zonally  aver¬ 
aged  model  to  changes  in  heating  is  severely  limited  by 
exclusion  of  baroclinic  instability  and  explicit  water  va¬ 
por  cycles.  Therefore  the  sensitivity  of  T  to  Wo  may  be 
primarily  a  model  tuning  issue.  However,  we  also  found 
that  T  is  sensitive  to  a  small,  directly  imposed,  strato¬ 
spheric  heat  source.  A  small  ~0.3  K  d-1  heat  source 
confined  to  the  tropical  lower  stratosphere  (Figure  lb) 
had  effects  on  both  T  and  w* ,  similar  to  those  obtained 
by  raising  the  tropospheric  heating  from  3  to  4  K  d_1. 
Variations  of  0.1-  0.3  K  d-1  in  stratospheric  net  heating 
such  as  obtained  in  experiment  7  (Figure  14)  are  well 
within  the  range  of  current  uncertainties.  It  has  been 
shown,  for  example,  that  tropical  lower  stratospheric 
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net  heating  rates  are  particularly  sensitive  to  details  in 
the  treatment  of  surface  albedo,  tropospheric  clouds, 
as  well  as  to  secondary  radiatively  active  species  such 
as  N2O,  NO2,  CH4,  and  O2  [Eluszkiewicz  et  al. ,  1996]. 
Furthermore,  the  stratospheric  aerosol  load  produced 
by  major  volcanic  eruptions  can  perturb  net  heating 
rates  by  0.2-0. 5  K  d-1  [ Kinne ,  1992;  Kinnison  et  a/., 
1994;  Eluszkiewicz  et  al.,  1996]. 

Thus  our  results  suggest  that  variability  in  strato¬ 
spheric  net  heating  may  be  a  significant  contributor  to 
observed  variations  in  upward  motion  in  the  tropics. 
Such  variability  may  have  to  be  considered  as  a  pos¬ 
sible  explanation  for  interannual  variations  in  the  rate 
of  vertical  transport  in  the  lower  stratosphere,  along 
with  the  variations  introduced  by  the  quasi-biennial  os¬ 
cillation  [Mote  et  a/.,  1996;  Randel  et  al. ,  1996].  A 
possible  connection  between  age  and  changes  in  dia- 
batic  heating  was  made  earlier  by  Schmidt  and  Khedim 
[1991]  to  explain  a  3  to  5  year  periodicity  in  values  of 
age  determined  from  balloon  borne  measurements  of 
CO2.  Variations  in  factors  such  as  cloud  cover,  sur¬ 
face  albedo,  or  actual  tropical  tropospheric  tempera¬ 
tures  may  be  related  to  El  Nino/Southern-Oscillation 
(ENSO).  Schmidt  and  Khedim  speculated  that  ENSO- 
related  effects  acting  on  the  meridional  circulation  could 
induce  the  3  to  5  year  periodicity  they  observe  in  age. 
Volcanic  aerosols  vary  in  an  unpredictable  manner,  but 
their  potential  impact  is  large.  Volcanic  aerosols  from 
Pinatubo  were  present  in  the  lower  stratosphere  dur¬ 
ing  much  of  the  early  UARS  observation  period  [e.g., 
Kinne  et  al. ,  1992;  Kinnison  et  al. ,  1994].  So  variations 
in  ascent  rates  may  be  present  in  the  UARS  HALOE 
and  MLS  data  as  result  of  this  changing  stratospheric 
aerosol  load. 

Appendix:  Idealized  Response  to  a 
Steady  Momentum  Source 

In  the  TEM  formulation  the  thermal  wind  relation 
is  used  with  (la)  and  (lb)  to  form  a  diagnostic  elliptic 
equation  for  the  residual  meridional  stream  function  ip* , 

Etw)  -  -(/  +  2S£) 

<A1> 

where  the  symbol  Ell()  denotes  a  complicated  elliptic 
operator,  which  depends  on  U  and  0  as  well  as  other 
zonal  mean  quantities.  The  term  in  (Al)  represents 
departures  from  geostrophic  balance,  which  are  param¬ 
eterized  as  discussed  by  Bacmeister  et  al.  [1995],  The 
component  stream  functions  ipx  and  and  associated 
velocities  discussed  in  section  4  are  simply  derived  from 
truncated  versions  of  (Al),  where  only  one  or  a  subset 
of  the  forcing  terms  on  the  RHS  has  been  retained,  for 
example, 


While  these  stream  function  components  contain  no  in¬ 
formation  that  is  not  already  present  in  their  source 
terms,  they  offer  a  straightforward  and  physically  mean¬ 
ingful  way  to  assess  the  relative  importance  of  the  nu¬ 
merous  drag  and  heating  mechanisms  in  a  2-D  model. 

In  this  Appendix  we  will  use  stream  function  com¬ 
ponents  to  examine  the  zonally  averaged  response  of  an 
idealized  atmosphere  to  a  localized  source  of  momentum 
[e.g.,  Haynes  et  al. ,  1991;  Holton  et  al.}  1995].  Haynes 
et  al.  showed  that  for  a  steady  momentum  source  the 
mean  meridional  mass  circulation  approaches  a  steady 
state  in  which  the  circulation  is  confined  along  lines  of 
constant  angular  momentum  below  the  forcing.  This 
configuration  of  momentum  forcing  and  residual  circu¬ 
lation  has  been  called  “downward  control.”  Here  we  will 
show  how  downward  control  is  established  by  a  reactive, 
thermally  driven  circulation,  which  arises  in  response  to 
temperature  changes  brought  about  by  the  initial  me¬ 
chanically  forced  circulation. 

We  specify  a  time  independent  momentum  forcing 
given  by 


The  peak  amplitude  of  the  steady  momentum  forcing 
A  Force  is  chosen  to  be  5  m  s'M"1  and  the  center  of 
the  forcing  is  chosen  to  be  at  2=45  km  and  0=45°  N. 
The  scales  for  the  forcing  are  set  to  A2Force=10  km  and 
A0Force=lO°.  The  net  heating  H  is  given  by  simple 
Newtonian  cooling  to  a  latitude  independent,  reference 
temperature  profile 

rj_^  _  0(0,  t)  —  ^Global (^)  (A3) 

'TNewt 

The  relaxation  time  r^ewt  is  constant  with  latitude  and 
height  and  is  chosen  to  be  4  days.  Parameters  in  (A2) 
and  (A3)  are  chosen  to  be  similar  to  those  used  by 
Haynes  et  al  [1991]. 

Figure  17  shows  the  evolution  of  the  meridional  mass 
stream  function  after  the  constant  momentum  forcing 
(A2)  is  turned  on  at  t= 0  in  an  atmosphere  initially  at 
rest.  The  top  row  of  panels  shows  the  time  evolution 
of  the  total  mass  stream  function  with  time  in¬ 
creasing  from  left  to  right.  The  middle  row  shows  the 
time  evolution  of  pxpx  (the  component  due  to  the  mo¬ 
mentum  forcing),  and  the  bottom  row  shows  the  time 
evolution  of  pxp^  (the  component  forced  by  thermal  re¬ 
laxation).  In  these  plots  a  maximum  in  pip*  implies 
clockwise  circulation.  Early  in  the  simulation  (after  6 
days)  the  total  mass  stream  function  pip*  is  close  to 
pip*x ,  implying  essentially  a  single  clockwise  cell  with  a 
broad  ascending  branch  on  the  equatorward  flank  and 
a  more  concentrated  descending  branch  on  the  pole- 
ward  flank.  The  initial,  thermally  forced  circulation 
pip^  is  weak.  However,  as  time  progresses,  the  temper¬ 
ature  field  is  displaced  farther  from  its  initial  equilib- 
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Figure  IT.  Residual  mass  stream  functions  pip  for  idealized  experiment  described  in  the 
Appendix.  Constant  body  force  centered  at  45° N  and  45  km  is  applied  starting  at  day  zero. 
Thermal  forcing  is  limited  to  Newtonian  cooling  to  reference  temperature  profile  with  a  time 
constant  of  4  days.  Top  row  shows  total  mass  stream  function  pipx-  Middle  row  shows  mass 
stream  function  forced  by  mechanical  forcing  alone,  pipx .  Bottom  row  shows  mass  stream  function 
forced  by  thermal  forcing  alone,  pipn.  Left-hand  column  shows  the  three  mass  stream  functions 
on  day  6,  middle  column  on  day  20,  and  right-hand  column  on  day  100. 


rium  state.  The  relaxational  thermodynamic  forcing  H 
becomes  stronger  as  does  the  thermally  forced  circula¬ 
tion  pip  ft .  The  thermally  forced  circulation  consists  of 
a  multiple  cell  pattern  with  a  strong,  horizontally  local¬ 
ized  clockwise  cell  below  the  forcing,  a  counterclockwise 
cell  developing  on  the  equatorward  flank  of  the  central 
cell,  and  a  third  weaker  counterclockwise  cell  developing 
on  the  poleward  flank. 

Eventually,  as  steady  state  is  achieved,  the  total  mass 
stream  function  pip*  approaches  the  “downward  con¬ 
trol”  configuration  discussed  by  Haynes  et  ai  [1991], 
with  motion  largely  confined  below  the  momentum  forc¬ 
ing.  This  can  be  seen  in  the  top  right-hand  panel  of  Fig¬ 
ure  17,  which  shows  pip *  after  100  days.  However,  the 
approach  to  the  steady  state,  downward  controlled  con¬ 
figuration  is  achieved  by  intensification  of  the  thermally 
forced  circulation  pip which  remains  highly  nonlocal. 
Horizontal  localization  (“downward  control”)  occurs  be¬ 
cause  pip  ft  almost  exactly  cancels  pipx  outside  of  the  re¬ 
gion  below  the  constant  momentum  forcing.  Through¬ 
out  the  simulation,  pipx  remains  approximately  con¬ 
stant  in  time  (Figure  17,  middle  row).  The  small  dif¬ 


ferences  in  pipx  between  day  6  and  day  100  are  due  to 
changes  in  the  coefficients  of  the  linear  elliptic  operator 
in  (Al)  which  result  as  U  and  9  evolve  with  time.  Nev¬ 
ertheless,  in  the  final  steady  state  both  pipx  and  pipft 
remain  highly  nonlocal  and  only  their  sum  exhibits  the 
horizontal  localization  required  of  the  state  steady  so¬ 
lution  by  the  downward  control  principle. 

Clearly,  in  the  simulation  just  described,  it  is  phys¬ 
ically  meaningful  to  speak  of  a  primary  residual  cir¬ 
culation  associated  with  momentum  forcing  and  a  sec¬ 
ondary  or  reactive  circulation  associated  with  thermo¬ 
dynamic  forcing/relaxation.  For  downward  control  to 
exist,  the  secondary  circulation  must  be  equal  and  op¬ 
posite  to  the  primary  circulation  in  regions  horizontally 
away  from  the  steady  momentum  source.  Thus  the 
eventual  magnitude  of  the  secondary  thermal  circula¬ 
tion  in  these  locations  is  completely  determined  by  the 
circulation  associated  with  the  momentum  forcing.  Im¬ 
mediately  below  the  imposed  momentum  forcing,  the 
situation  is  more  complicated,  P^n  becomes  more  in¬ 
tense  than  pnp*x .  Thus  the  total  stream  function  below 
the  imposed  forcing  winds  up  closer  in  magnitude  and 
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shape  to  the  thermally  forced  stream  function,  despite 
the  fact  that  its  ultimate  “cause”  is  the  mechanically 
forced  circulation  in  the  middle  row  of  Figure  17. 

An  implication  of  the  simulation  in  Figure  17  is  that 
for  relaxational  net  heating  the  adiabatic  response  to  a 
distant  momentum  source  represents  a  reasonable  esti¬ 
mate  of  the  maximum  impact  that  momentum  source 
can  have  on  the  local  circulation.  In  other  words,  the 
thermally  reactive  circulation  (Figure  17,  bottom  row) 
will  be  at  most  of  a  similar  magnitude  to  the  adiabatic 
response  (Figure  17,  middle  row)  horizontally  well  away 
from  the  momentum  forcing.  Thus  it  seems  reasonable 
to  use  the  adiabatic  stream  function  response  to  a  given 
momentum  forcing  to  assess  the  magnitude  of  its  pos¬ 
sible  remote  effects.  We  would  argue  therefore  that  the 
small  contribution  of  to  the  total  vertical  velocity 
in  the  tropical  stratosphere  (Figure  12),  argues  against 
a  significant  impact  by  midlatitude  drag. 

Information  obtained  by  decomposing  the  stream  func¬ 
tion  into  components  does  not  contradict  results  ob¬ 
tained  by  steady  state  reasoning  where  steady  state 
conditions  apply.  In  the  example  just  examined,  the 
final  net  heating  away  from  the  imposed  momentum 
source  approaches  zero,  so  that  w*  =  0  there  is  con¬ 
sistent  with  the  steady  state  solution.  However,  the 
zero-motion  state  away  from  the  momentum  forcing  is 
achieved  by  two  nonzero  circulations  associated  with 
the  thermal  forcing  and  the  momentum  forcing  of  the 
elliptic  equation  (Al).  Examining  component  stream- 
fucntions  may  be  helpful  in  situations  where  forcings 
vary  with  time  and  relaxational  time  constants  are  dif¬ 
ficult  to  determine.  Furthermore,  even  in  steady  state 
situations,  component  stream  functions  may  be  helpful 
in  evaulating  the  sensitivity  of  circulations  to  changes 
in  heating  and  drag  parameters. 
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